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ABSTRACT 

This  study  employs  a  radiational  model,  with  large-scale  cloud 
parameterization  in  several  layers,  in  order  to  compute  the  absorption 
of  solar  insolation  by  the  ocean's  surface  and  by  atmospheric  layers, 
and  to  compute  the  planetary  albedo  for  the  FNWC  primitive  equation 
model.   Solar  insolation  dispositions  were  computed  on  a  gridpoint 
basis  using  water-vapor  mass  and  cloud  amounts.   Using  emissivity  for- 
mulas after  Sasamori,  the  long-wave  cooling  effects  were  calculated  at 
the  earth's  surface  and  over  the  tropospheric  layers,  and  were  found 
to  be  dependent  upon  the  cloud  parameterization. 

FNWC  data  over  oceanic  gridpoints  for  16  July  1974  were  used  to 
test  two  forms  of  the  cloud  parameterization.   The  objective  was  to 
determine  the  parameterization  which  better  verified  the  radiational 
balance  at  the  top  of  the  atmosphere  as  a  function  of  latitude  as 
compared  with  satellite  climatology.   The  smaller  cloud  parameteriza- 
tion gave  the  better  verification. 
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I.   INTRODUCTION 

This  thesis  is  a  study  of  a  radiative  heating  parameterization  for 
use  in  the  Fleet  Numerical  Weather  Central  (FNWC)  prediction  system. 
The  study  has  as  a  primary  objective,  the  analysis  of  the  radiational 
and  heat  balance  of  the  ocean-atmosphere  system  utilizing  FNWC  gridded 
data-fields  at  constant  pressure  levels  for  16  July  1974.   The  grid- 
point  data  were  selected  along  four  oceanic  meridians,  three  in  the 
Pacific  and  one  in  the  Atlantic,  with  most  of  the  gridpoints  located 
in  the  Northern  Hemisphere  (Fig.  1) . 

The  specification  of  amounts  of  clouds  in  two  designated  layers 
has  the  most  influence  on  the  radiative-model  dispositions  (short-  and 
long-wave) .   The  initial  specification  of  the  fractional  amounts  of 
CL(1)  and  CL(2)  were  based  on  large-scale  formulations  developed  by 
Smagorinsky  (1960)  and  used  in  a  similar  study  based  on  the  data  day 
16  January  1974  by  Spaeth  (1975) . 

The  governing  equations  of  the  radiational  transfer  of  the  model 
have  been  derived  by  Martin  (1972,  1974) ,  who  modified  both  the  solar 
and  terrestrial  radiation  transfers  to  respond  to  the  presence  of 
clouds.   These  clouds  were  of  specified  amounts  in  the  two  layers,  one 
of  which  is  a  mid-level  cloud  and  the  other  a  low-level  cloud.   The 
radiational  model  in  use  here  has  similarities  to  those  in  use  in  UCLA 
and  NCAR  General  Circulation  models. 

The  physical  description  of  the  radiative  model  may  be  applied  to 
any  scale  of  analysis  for  which  there  is  adequate  resolution  of  the 
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temperature  and  moisture  data  in  the  vertical.   In  the  horizontal,  the 
reliability  of  the  data  used  here  is  consistent  with  that  of  the  FNWC 
analysis  to  gridpoints,  and  is  typically  reported  to  the  nearest  tenth 
of  a  degree  with  regard  to  temperature  and  dewpoint.   The  radiational 
computations  made  here  are  applied  to  FNWC  gridpoints  and  are  designed 
to  make  a  one-hour  forward-time  step  applicable  to  the  FNWC  primitive 
equation  forecast  model,  with  special  adaptations  to  their  a-levels. 

The  Smagorinsky  (1960)  cloud-specification  of  CL(1)  and  CL(2)  when 
used  with  the  model  of  this  study  resulted  in  global  albedo  estimates 
that  were  too  high  especially  in  the  tropics  as  compared  to  the  satel- 
lite data  of  Raschke  (1973)  ,  thus  leading  to  under-estimates  of  the  net 
radiative  balance  at  the  top  of  the  atmosphere  and  at  the  ocean  surface, 
This  could  be  a  result  of  inadequately  resolved  relative  humidities, 
which  in  reality  might  well  be  lower  than  inferred  by  the  FNWC  tropical 
analysis. 

The  comparative  results  afforded  by  the  2/3-CL  parameterization 
gave  reasonably  close  agreement  with  the  radiative  climatology  of 
Raschke  et  al  (1973),  for  the  NIMBUS  III  period  16-31  July  1969. 
Oceanic  and  tropospheric  balances  were  computed  for  both  sets  of  cloud 
conditions  (full-CL  and  2/3-CL)  for  16  July  1974.   It  was  also  neces- 
sary to  augment  the  purely  radiative  model  by  including  a  turbulent 
boundary-layer  model  for  sensible  and  latent  heat  transports  over  the 
ocean  as  adapted  from  Kaitala  (1974) . 

The  radiation  package  results  computed  with  the  2/3-CL  cloud  model 
compared  very  favorably  with  satellite  climatology,  though  only  one 
data  day  (16  July  1974)  was  used  as  representative  of  mid-July  sound- 
ings over  the  ocean,  both  for  the  model  and  the  Raschke  satellite 
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climatology.   However,  the  latter  climatology  was  based  on  the  data 
period  16-31  July  1969. 

The  major  effect  noted  by  the  use  of  the  2/3-CL  as  contrasted  with 
full-CL  parameterization  was  the  reduction  in  global  albedo  and  the 
resulting  greater  contribution  to  the  surface  net  heating  rates.   This 
result  is  in  general  agreement  with  recent  observations  by  Von  der  Haar 
and  Hanson  (1969) ,  Von  der  Haar  and  Oort  (1973)  and  Raschke  et  al 
(1973) .   However,  the  major  discrepancies  in  the  radiative  model  re- 
sulted from  the  specification  of  layered  cloud  amounts  in  the  tropics 
(in  a  manner  identical  to  that  used  poleward  of  the  tropics)  where 
the  model  gave  somewhat  excessive  albedos.   It  is  felt  that  the  geo- 
metry of  solar  radiation  streams  impinging  upon  subgrid  sized  cumulus 
cells  in  the  tropics  should  be  remodeled  so  as  to  divert  a  greater  per- 
centage of  solar  radiation  downward,  as  compared  with  the  characteris- 
tic cloud-layering  and  resultant  reflective  effects  which  seem 
permissible  in  middle  and  high  latitudes. 
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II.   DATA  PREPARATION 

A.   INITIAL  DATA  FIELDS 

The  temperature  and  humidity  data  used  in  this  study  were  arranged 
in  the  form  of  soundings  taken  along  four  oceanic  meridians  (Fig.  1) 
of  the  Fleet  Numerical  Weather  Central  (FNWC)  Northern  Hemisphere  ana- 
lyses on  16  July  1974.   Computation  of  radiational  dispositions  and 
of  other  heat  budget  terms  were  made  at  these  gridpoints.   Oceanic 

locations  for  these  computations  were  chosen  because: 

p 

1.  The  constant  a  surfaces  (where  0  -  —  )    of  the  FNWC  primitive 

TT 

equation  system  are  close  to  being  constant  pressure  levels. 

2.  The  maritime-area  heating  rate  computations  are  likely  to  be 
representative  of  the  month  of  July  1974  as  contrasted  with  computa- 
tions made  for  16  July  1974  for  a  corresponding  set  of  gridpoints 
over  land. 

The  three  meridians  (and  their  respective  number  of  soundings) 
selected  over  the  Pacific  Ocean  were  located  at  125W  (25  soundings) , 
170W  (25  soundings)  and  145E  (17  soundings) .   The  Atlantic  Ocean 
meridian  was  35W  (26  soundings) .   This  method  of  selecting  "sound- 
ings" along  the  indicated  meridians  of  the  FNWC  polar  stereographic 
map,  made  it  unnecessary  to  employ  spatial  interpolation  between 
original  data  gridpoints  along  the  meridians.   Data  along  line  3  in 
the  Pacific  was  not  extended  southward  of  gridpoint  (9,55)  because 
they  fell  over  land  masses  (New  Guinea  and  Northern  Australia)  where 
the  surface  temperatures  and  other  sounding  features  were  unrepre- 
sentative of  the  oceanic  values.   For  similar  reason  two  other 
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|1.43) 


(t.32) 


(63.1) 


Figure  1.   FNWC  polar  stereographic  grid  and  meridians  (lines  1,  2, 
3,  and  4)  selected  for  study.   The  longitudes  A  are  shown 
for  each  meridian  as  well  as  the  extent  considered  of  each 
meridian. 
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gridpoints  at  (24,24)and  (56,6)  which  also  fell  just  over  land  masses, 
were  replaced  with  more  representative  oceanic  soundings  from  the 
gridpoints  (23,24)and  (56,7) ,  respectively. 

The  gridpoint  soundings  in  the  form  of  Table  1(a)  were  taken  from 
the  original  FNWC  63-by-63  surface  analysis,  as  well  as  the  vertical 
distribution  of  T(P)  at  ten  additional  levels  ranging  up  to  100  milli- 
bars.  Moisture  parameters  were  provided  as  the  vapor  pressure  at  the 
top  of  the  constant  flux  level,  and  as  the  dewpoint  depression  at  five 
selected  levels  between  the  925  millibar  and  400  millibar  levels. 
These  original  soundings  were  then  modeled  into  a  radiative  sounding 
(Table  I (b) )  corresponding  in  the  vertical  to  FNWC  primitive  equation 
prediction  levels  up  to  100  mb  (see  Fig.  2) . 

Data  from  0000GMT,  16  July  1974  were  used  for  the  Pacific  sound- 
ings while  1200GMT,  16  July  1974  data  were  used  for  the  Atlantic  sound- 
ings.  These  times  were  used  so  that  the  analysis  times  would  correspond 
as  closely  as  possible  to  local  solar  noon  in  each  of  the  indicated 
areas. 

Since  the  standard  instrument  level  vapor  pressure  (e  .  )  was  miss- 
air 

ing  from  the  original  FNWC  analysis  at  all  gridpoints,  e  .   was  approxi- 

ai  r 

mated  using  the  FNWC  field  of  e  ,  a  computed  value  of  vapor  pressure 

at  about  20  meters  above  MSL  in  the  turbulent  boundary  layer0   FNWC 

values  of  T   and  e   over  the  ocean  had  been  made  available  using  an 
xx 

operational  planetary  boundary  layer  model  detailed  by  Kaitala  (1974) . 
To  perform  radiative  calculations,  it  is  necessary  to  have  water- 
vapor  and  CO   absorber  masses  and  cloud  amounts,  CL(1)  and  CL(2) 
within  certain    k-level  boundaries  (Fig.  2).   All  soundings  in  this 
study  start  at  sea  level  with  the  approximation  of  surface  pressure 
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Table  1(a).   Example  of  typical  FNWC  sounding  for  gridpoint  (2,2). 
The  humidity  parameters  between  925,  ...,400  mb  are  dewpoint  depres- 
sions, whereas  the  values  at  the  surface  and  the  top  of  the  constant 
flux  layer  (9999*)  are  vapor  pressures. 


Pressure 

(mb) 

T(°C) 

Surface, 

1000  mb 

26.1 

925 

17.5 

850 

13.7 

700 

4.8 

500 

-11.8 

400 

-23.4 

300 

-38.2 

250 

-47.1 

200 

-56o9 

150 

-67.8 

100 

-81.0 

9999* 

T   =  23.5 

X 

Humidity  Parameters 


e 

X 

=  23.8  mb 

5.4  °C 

7.7 

12.7 

13.5 

15.3 

*  Code  9999  indicates  data  taken  from  the  top  of  the  constant  flux 
level  (CFL)  of  the  FNWC  initial  data  program. 
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Table  I  (b) .   Example  of  the  corresponding  radiative  sounding  with  mix- 
ing ratio  listed  at  odd  k-levels  (Fig.  2).   Additionally,  water-vapor 
and  CO   absorber  masses,  and  cloud  amounts  CL(1)  and  CL(2)  (after 
Smagonnsky  (I960))  are  also  listed,  as  these  parameters  have  been 
modeled  in  the  radiative  theory  presented  in  this  study. 


Mixing 

Absorber  M 

Pressure 

Ratio 

Water  Vapor 

(mb) 

T(°C) 

g/kg 

(gm/cm^) 

1000 

26.1 

14.93 

900 

15.7 

8.45 

800 

10.9 

0.0 

1.58 

700 

4.8 

3o04 

600 

-2.4 

0.0 

2.06 

500 

-11.8 

1.02 

400 

-23.4 

0.36 

2.18 

300 

-38.2 

0.15 

200 

-56.9 

0.0 

2.20 

100 

-81.0 

0.00 

2.20 

0 

-81.0 

0.0 

2.20 

C09 
(cm/ cm  ) 


45.53 

83.53 

113.35 

133.99 
138.67 
143.35 


Smagorinsky 
Cloud  Amounts 
CL(1)  ,  CL(2) 


0.24,  CL(2) 


0.00,  CL(1) 
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T. 


Layer  (0,2) 


(2,4) 


(4,6) 


(6  8) 


J, 

(8,10) 

1 


Level- 
Idenlifers 


U(O),  C(0) 


O-level 
O  =  0 


T,q 0.1 


2- 


U(2),  C(2) 


0=0.2 


0.3 


0=0.4 


o=0.6 


7 T,q 0.7 


U(8),  C(8) 


0=0.8 


10 


U(10),  C(10) 


o=1.0 


Figure  2.   Five-layer  radiative  sounding  used  in  this  study.   Levels 
are  identified  by  their  values  on  the  k-scale,  while  layers  are  identi- 
fied by  their  level  boundary  indices  in  parentheses,  e.g.  (8,10). 
Pressure-scaled  water  vapor  and  CO  mass  increments  Au  and  AC,  respec- 
tively are  integrated  relative  to  the  surface  and  the  resulting  U  and  C 
are  carried  at  even  levels.   The  temperature  T  is  retained  at  all  levels, 
Amounts  of  clouds  CL(1)  and  CL(2)  in  the  layers  shown  have  been  para- 
meterized for  consideration  of  their  radiative  effects. 
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7r  =  1000  mb.   Therefore,  the  k-levels  are  associated  approximately  with 
the  FNWC  levels  P   =  1000.,  900.,  800.,  ....,  200.,  100.,  0.0  mb  respec- 
tive  too  -  1.0,  0.9,  ...,  0.1,  0.0. 

B.   INTERPOLATIVE  PROCESSING  TO  K-LEVELS  IN  RADIATIVE  SOUNDINGS 

1.  Temperature  Profiles 

The  gridpoint  temperatures  were  listed  at  each  mandatory  level 
of  Table  1(a)  between  1000.,  ...,  100  mb.   The  temperature  was  assumed  to 
be  isothermal  from  100  mb  to  p=0.0  mb.   The  temperature  T   was  set  equal 
to  the  FNWC  listed  sea-surface  temperature.   The  radiative  sounding  tem- 
peratures for  the  remaining  k-levels  were  obtained  from  either  their 
corresponding  listed  temperature  value  or  by  a  three-point  Lagrangian 
interpolation  scheme  (Eq.  2-1,  (Spaeth,  1975)),  to  level  k  when  the 
listed  temperature-profile  does  not  include  the  value  T  . 

2.  Moisture  Profiles 


The  moisture  from  the  original  FNWC  soundings  (Table  I  (a) )  were 

converted  into  mixing  ratios  at  each  of  the  original  sounding  levels. 

The  near-surface  vapor  pressure  (taken  as  e  )  was  used  to  calculate  the 

x 

mixing  ratio  at  k=10,  as  described  by  Spaeth  (1975) : 


qn^  =  621.97  (e  /1000)  (2-1) 

10  x 


The  remaining  mixing  ratios  q  were  calculated  from  the  original  FNWC 
sounding  (Table  1(a)),  using  a  dewpoint  depression  formula  (Spaeth,  1975, 
(Eq.  2-3)).  The  computed  q-values  were  subsequently  interpolated  to  k- 
levels  using  the  previously  mentioned  three-point  Lagrangian  procedure 
applied  to  successive  q-values  in  the  original  sounding.   Resulting 
values  of  the  interpolated  q-values  are  shown  for  the  case  of  the  radia- 
tive sounding  at  gridpoint  (2,2)  (Table  1(b)). 
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Due  to  the  fact  that  most  radiosonde  humidity  data  for  p<300  mb 
is  either  unreliable  or  not  available,  a  power- law  extrapolation  for- 
mula (Spaeth,  1975) 

Vq500=    <V50°)X  (2"2) 

was  used  to  obtain  q-values  at  k=3f2,l. 

This  extrapolative  procedure  for  obtaining  q-values  at  high 
levels  was  first  suggested  by  Smith  (1966) ,  utilizing  the  least  squares 
estimation  of  X   according  to 

6 


(2-3) 


with  y.  =  log  (q./qc)  and  x.  =  log  (p  . /p  )  .   The  Smith  method  was  then 

further  tested  by  Spaeth  (1975)  for  the  vertical  scale  of  the  FNWC 

initial  data  soundings.   The  resulting  correlation  coefficient  R 

yx 

was  found  here   to  lie  in  the  range  .95  to  .99.   Such  high  correlations 
indicate  that  the  procedure  described  by  Spaeth  (1975)  for  determining 
upper  atmosphere  q-values  has  realistic  estimation  value  for  the  verti- 
cal scale  involved  in  the  radiative  sounding 0 
3.   Pressure-Scaled  Absorber  Masses 

The  pressure-scaled  water  vapor  absorber  mass  Au  in  a  layer 
(Fig.  2)  were  calculated  for  the  five  odd  k-levels  by  using  the  computed 
mixing  ratio  values  (Eq.  2-7,  Spaeth,  (1975)).   The  equation  for  com- 
puting the  integrated  pressure-scaled  water  vapor  mass  is  given  by  the 
algorithm  described  by  Eq.  2-8  of  Spaeth. 
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A  similar  routine  was  used  in  the  computation  of  the  carbon- 
dioxide  scaled  absorber  mass  (Spaeth,  1975,  Eqs .  2-10,11).   To  be  noted 
here  is  the  fact  that  CO  absorber  masses  have  been  stated  in  terms  of 
N.T.P.  "pressure-scaled"cms,  since  the  CO   absorption  coefficients  are 
generally  stated  in  terms  of  the  reciprocal  of  this  CO  mass  unit. 

C.   CLOUD  PARAMETERIZATION 

The  relative  humidities  and  thus  the  saturation  vapor  pressure  at 
levels  k=5  and  k=9  are  used  in  the  calculations  of  the  fractional  cloud 
cover  CL  and  CL   in  layers  (4,6)  and  (8,9)  respectively,  as  depicted 
in  Fig.  2.   From  the  computed  relative  humidity  values  at  k=5  and  k=9, 
the  fractional  cloud  amounts  for  each  layer  were  parameterized  using 
the  equation  (after  Smagorinsky,  (I960)): 

CL(1)  =  2.0  (RH(5))  -  0.7  (2-4 (a) ) 

CL(2)  -  3.33(RH(9))  -  2.0  (2-4 (b)) 

Smagorinsky ' s  parameterization  of  CL(1)  and  CL(2)  permits  cloudiness 
fractions  of  1.30  and  1.333  respectively  with  RH=1.0.   Values  of  CL 
greater  than  1.0  were  considered  initially  by  Smagorinsky  (1960)  to 
suggest  the  amount  of  supersaturation  which  accompanies  precipitation. 

All  radiative  models  in  use  at  the  present  time  limit  CL  in  accord- 
ance with  0  <  CL  <  1.0.   In  the  present  radiative  package  there  is  little 
opportunity  to  deduce  precipitation  and/or  supersaturation  amounts. 
Hence,  the  formulations  of  CL  and  CL  of  Eqs.  2-4  were  limited  to  a 
maximum  of  1.0  in  each  layer  in  the  "full-CL"  case.   In  addition,  the 
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Smagorinsky  computations  for  CL  and  CL  were  reduced  to 


CL  »  =  2/3  (CI^)  (2-5  (a)) 

CL2'  =  2/3  (CL2)  (2-5 (b) ) 


regardless  of  whether  the  original  CL  values  exceeded  1.0  or  not,  so 
that  the  revised  cloud-cover  amounts  could  be  considered  for  purposes 
of  estimating  shifts  in  the  global  radiative  balance.   Computations 
made  with  cloud  parameterization  of  Eq.  2-5,  are  termed  the  2/3-CL  case. 

The  1/3  reduction  inherent  in  CL  * ,  CL  '  relative  to  Smagorinsky ' s 
CL  ,  CL  were  initiated  in  this  study  in  an  attempt  to  tune  cloud 
amounts  to  give  radiation  results  in  closer  agreement  with  recent  satel- 
lite climatology  of  Raschke  et  al,  (1973).   This  satellite  climatology 
suggests  the  use  of  smaller  cloud  amounts  than  that  specified  by 
Smagorinsky.   Fractional  cloud  amounts  were  considered  to  be  functions 
of  large  scale  effects  only.   Therefore,  small  scale  convective  acti- 
vity, seasonal  conditions  (except  for  the  soundings)  and  latitudinal 
effects  were  not  considered  in  specifying  the  reduction  factor  of 
Eq.  2-5.   The  reduced  cloud-parameterization  was  introduced  here  for 
estimating  large  scale  radiational  effects  only. 

D.   CLOUD-AREA  COVERAGES 

Since  Eqs .  2-4a,b  or  2-5a,b  gave  the  fractional-coverage  of  the 
gridpoint  area  by  the  appropriate  cloud- type,  the  gridpoint  area  may 
be  thought  of  as  broken  into  random  fractional  segments  of  size 


W(0,0)  =  (1-CL1)*(1-CL  )  (2-6(a)) 
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wherein  there  is  a  combination  of  clear-over-clear  segments  in  the 
layers.   Similarly,  the  gridpoint  area  has  the  fractional  area  of  cloud 
coverage 


W(l,l)  =  CL   *  CL2  (2-6 (b)) 


of  upper-cloud  amount  CL  overlying  lower-cloud  amount  CL  .   Likewise 
the  combinations  of  cloudy-over-clear  and  clear-over-cloudy,  by  layers, 
may  be  visualized  as  occurring  with  the  weights 


W(1,0)  =  CL±    *  d-CL2)  (2-6(c)) 


and 


W(0,1)  =  (1-CL1)  *  CL2  (2-6(d)) 


Regardless  of  whether  the  full-cloud  amounts  CL  ,  CL  of  Smagorin- 


sky  (Eq.  2-4) ,  or  the  2/3-CL  amounts  of  Eq.  2-5  were  utilized,  it  was 
useful  to  carry,  for  radiation  computations  of  each  sounding,  the 
relative  weights  or  fractions  of  the  gridpoint  area  exposed  to  the 
specified  cloud- layer  combinations.   Henceforth,  the  symbols  are  denoted 
by  W(0,0),  W(l,l),  W(1,0),  W(0,1)  and  as  given  by  Eqs .  2-6(a,b,c,d) 
suggest  overcast  (1)  or  clear  (0)  cloud  amounts  in  the  indicated  layers 
(Fig.  2) ,  the  first  index  referring  to  the  level  of  CL   and  the  second 
to  that  of  CL  .   The  usefulness  of  this  computational  device  will  be 
exemplified  in  Sections  III  and  IV,  where  the  procedures  for  the  terres- 
trial and  solar  radiational  computations  are  discussed  and  the  results 
are  summarized  over  the  full  set  of  soundings. 
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III.   TERRESTRIAL  RADIATION 

A.   THEORETICAL  AND  EMPIRICAL  BASIS 

Empirical  formulas  were  developed  by  Sasamori  (1968)  for  flux  emis- 
sivities  in  the  atmosphere  as  associated  with  computations  required  for 
the  heat  balance  requirement  of  the  NCAR  General  Circulation  Model. 
Sasamori  derived  the  empirical  emissivity  formulas  for  water  vapor  and 
CO  by  comparison  with  the  theoretical  values  built  into  the  Yamamoto 
Radiation  Chart  (1952) .   This  chart  has  proved  to  be  quite  accurate  for 
computational  checking  of  the  Sasamori  emissivities  and  was  used  by 
Spaeth  (1975)  and  Warner  (1974)  as  a  systematic  guide  for  integration 
of  the  radiative  transfer  formula  developed  by  Martin  (1972,  1974),  who 
adapted  the  Sasamori  formulas  to  the  particular  layers  of  interest  in 
gridpoint  computations  of  the  FNWC  primitive  equation  model  (Fig.  2). 

The  essential  long-wave  flux  formulas  required  for  use  in  the  FNWC 

heating  package  are  the  following: 

* 

F    =  net  IR  flux  at  earth,  k=10 

* 
F.    =  net  IR  flux  at  level  k=6 
6 

* 

F    =  net  IR  flux  at  level  k=2 

F610  =  net  IR-flux  divergence  in  the  layer  (6,10) 
F26   =  net  IR-flux  divergence  in  the  layer  (2,6) 

To  get  the  IR-flux  divergences  in  the  layers  (6,10)  and  (2,6),  the 

*      *       *      * 
differences  F.   -  Fn^   and  F„   -  F^   must  be  computed.   The  detailed 
6      10        2      6 

scheme  for  making  such  computations  with  various  combinations  of  cloud 
cover  CL(1)  and  CL(2)  is  similar  to  that  developed  by  Martin  (1974)  and 
as  reproduced  in  detail  by  Spaeth  (1975) . 
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In  order  to  make  IR  net-flux  calculations  along  the  path  of  integra- 
tion, there  must  be  a  physically  sound  representation  of  the  emissivity 

(e   )  as  a  function  of  both  water-vapor  and  CO  absorber  masses  in  layers 
wc  2 

along  the  sounding.   For  a  further  discussion  of  the  emissivity  formulas 
used  in  the  quadrature  scheme,  refer  to  Spaeth's  Appendix  A  (1975). 

B.   NET  FLUX  FORMULATION 
1.   At  Level  k=10 

The  radiative  sounding  as  depicted  in  Table  1(b),  was  computed 
as  the  combination  of  parameters  U(k,10),  C(k,10)  and  T   for  each  re- 
quired  level,  k=10,8 , . . . , 2 , 1 ,0.   Cloud  parameters  CL(1)  and  CL(2)  are 
also  listed  at  each  gridpoint  and  in  general  are  non-zero.   The  grid 
area  was  then  considered  to  be  composed  of  areal  fractions  (weights) 
defined  in  Eqs .  2-6(a,...,d)  and  denoted  by  the  symbols  W(0,0),  W(l,l), 

W(1,0)  ,  W(0,1)  . 

* 
The  composite  net  flux  F     (CL  ,  CL  )  at  level  k=10  at  each 

gridpoint  is  then  constructed  by  using  the  appropriate  weight  factor  to 

* 
multiply  the  corresponding  reference  net  flux  F    computations , defined 

*  *  * 

for  the  special  cloud-cover  cases, F    (0,0),  F    (1,0),  F    (0,1), 

* 

F    (1,1)  .  As  a  result,  the  composite  net  flux  may  be  written 


F10  (CL1'CL2)  =  W(°'°)F10  (0,0)+W(l,0)F10*(l,0)+W(0,l)F10*(0,l) 

+W(1,1)F10*(1,1)  (3-1) 


Spaeth  (1975)  has  listed  these  reference  net  flux  formulations  in  his 
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Eqs.  3-6,7,8.   Using  the  definitions  of  W(0,0),  W(1,0),  W(0,1),  W(l,l) 

* 
F    (CL  ,CL  )  can  be  shown  to  assume  the  form 


F10  (CLrCL2)  =  [1-CL(2)]{(B10-B6)-.5[e:wc(8,10)(B10-B8) 

+  (£   (8,10)+£   (6,10)) (B-B.)]} 
wc        wc         8   6 

+  (1-CL(2)) (1-CL(1)){B  -.5[(£   (6,10) 

6      wc 

+  £   (4,10))  (B  -B.)+(£   (4,10)  (3-2) 

wc         6   4    wc 

+  £   (2,10))  (B  -B  )+(£   (2,10) 
wc  4   2     wc 

+  £   (1,0))  (B  -BJ+E   «0(10),T  )*BJ} 
wc        2   1   wc        1    1 

+  CL(2){(B   -BQ)  [1-.5E   (9,10)]}  . 
10   9       wc 


Here 


B,  =  1.170403  x  10  ?  T,  4  (3-3) 

k  k 


is  the  Stefan-Boltzmann  blackbody  flux  in  ly  day  '"  at  temperature  T  . 

K. 

Further,  £   (U,  ,C,  ,  10)  is  the  combined  water-vapor  and  CC< 
wc   k   k  2 

emissivity  along  the  path  from  level  10  to  level  k.   This  emissivity  is 
considered  by  Sasamori  to  be  temperature  independent  for  T  >  210K, 
whereas  £   represents  the  temperature  dependent  emissivity  applicable 
for  T  <  210K  (see  pp.  136-137,  Spaeth,  1975). 

The  reference  net  fluxes  F    of  Eq.  3-1  are  associated  with  (1) 
clear  skies  in  both  layers,  (2)  overcast  in  the  upper  layer  only,  (3) 
overcast  in  the  lower  layer  only  and  (4)  overcast  in  both  layers, 

respectively. 

* 

2.   Net  Flux  F 

6 — 

* 

The  formula  for  F   (CL. ,CL  )  has  been  developed  by  Spaeth  (Eq. 
6     12 

3-10,  197  5)  in  a  manner  analogous  to  the  derivation  of  the  weighted 
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* 
F    .   The  result  is  reproduced  for  completeness  after  (pp.  48,  Spaeth, 

1975). 

F6*  -  [l-CL(l)]{B8-.5[Cwc(6,8)(B8-B6)+E:wc(4,6)(B6-B4) 

+  (£wc(4'6)+ewc(2'6))  (VB2)+(ewc(2'6)+Cwc(1'6))  (B2'B1) 

+  6   ((0,6),T1)*Bl}+(l-CL(l))  (l-CL(2)){(B1  -BQ)  [1. 
WC        1    1  10   8 

-.5(£   (6,8)+E   (6,10))]}+CL(1){(B  -Bj* 
wc       wc  8   6 

[1-.5E   (6,8)] }+CL(l) (1-CL(2)){(B   -BQ)* 
wc  10   8 

[1-.5(E   (6,8)+e   (6,10))]}  .  (3-4) 

wc       wc 


* 

3.   Net  Flux  F 

The  net  flux  at  level  k=2  can  be  calculated  in  a  similar  fashion 

*       * 

to  F,.   and  F.    ,  with  the  formulation  of  the  individual  reference  net 
10       6 

*  * 

fluxes  F    as  previously  described.   The  formula  for  F   (CL  ,CL  )  is 

*  * 

developed  analogous  to  that  for  F    (CL  ,CL  )  and  F   (CL  ,CL  ) . 

F  *  =  [1-CL(1)]{B  -,5[£   (2,4)  (B  -B_)  +  (e   (2,4) 
2  8     wc        4   2     wc 

+   £       (2,6))  (B-BJ  +  (£   (2,6)+£   (2,8))* 
wc        6   4    wc       wc 

(BQ-Bc)+£   (1,2) (B _-B.)+E   ((0,2),  T.)*B.)]} 
8   6   wc       2   1   wc         11 

+  (l-CL(l))  (1-CL(2))*{(B.  -BQ)  [1-.5(E   (2,8) 

10   8        wc 

+  £   (2,10))]}+CL(1){B  -.5[C   (2,4)(B-B„) 
wc  4     wc       4   2 

+£   (1,2) (B-B.J+E   ((0,2),  T.)*B .]}  (3-5) 

wc       2   1   wc         11 

A  new  parameter,  the  total  outgoing  long-wave  radiation  to  space 

(FF2)  can  be  readily  defined  from  the  expression  for  F   by  setting  to 

zero  the  terms  of  Eq.  3-5  representing  the  downward  flux  through  level 

k=2.   These  terms  are  just  those  denoted  by  the  symbols  e   ( (0, 2) ,Tn ) *Bn 

wc        1   1 
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and  e   (1,2)*(B  -Bn ) .   If  the  long-wave  radiative  model   including  its 
wc         2   1 

cloud  parameterization   is  realistic,  the  modeled  outgoing  flux  to  space 
(FF2)  should  compare  reasonably  with  that  observed  by  NIMBUS  III  satel- 
lite during  the  same  period,  after  Raschke  et  al  (1973) .   The  most 
nearly  comparable  satellite  climatological  period  was  16-31  July  1969. 
These  satellite  observations  were  obtained  from  the  NIMBUS  III  satellite 
atlas  of  Raschke  et  al  (1973)  for  the  July  dates  indicated  and  used  for 
comparison  with  the  model  values  of  FF2,  point  for  point  (at  5°  latitude 
intervals)  along  the  same  four  meridians  of  Fig.  1.   Also  the  latitu- 
dinally-distributed  mean  satellite  fluxes  were  compared  with  the  corres- 
ponding averaged  values  of  FF2.   The  results  showing  the  comparisons  of 
the  latitudinally-averaged  distribution  of  the  model  computations  of 
FF2  and  the  corresponding  satellite  observations  are  listed  by  latitude 
in  Table  V. 

C.   APPLICATIONS  TO  HEAT  BALANCE  COMPUTATIONS 
1.   At  Upper  and  Lower  Boundaries 

In  order  to  compute  the  heat  balance  at  the  top  of  the  earth- 
tropospheric  system  the  composite  F   and  the  total  insolation  absorbed 
below  level  k=2  at  each  gridpoint  are  required.   Moreover,  a  radiative 

balance  is  computed  at  the  surface  for  each  gridpoint  using  the  earth's 

* 
absorbed  insolation  (Section  IV)  and  the  composite  net  flux  loss  F    . 

In  order  to  consider  the  possibility  of  a  heat  balance  at  the 

earth's  surface,  rather  than  merely  a  radiative  balance,  it  is  necessary 

to  add  a  term  representing  combined  sensible  and  latent  heat  transfer 

across  the  air-sea  interface  (Section  V) .   This  latter  transfer  rate  at 

each  gridpoint  was  adapted  from  the  FNWC  primitive  equation  model,  after 
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the  discussion  of  Kaitala  (1974)  and  based  upon  the  FNWC  data  of  16 
July  1974. 

2.   Intermediate  Levels 

Consideration  of  the  radiation  balance  in  the  atmosphere  re- 
quires computation  of  the  long-wave  cooling  effects  caused  by  the 
flux  divergences  F26  and  F610.   F26  and  F610  are  defined  as 

F26  =  F2*  -  F6*  (3-6) 

F610  =  F6*  -  F1Q*  (3-7) 

where  the  symbols  (2,6)  and  (6,10)  indicate  the  layer  boundaries  in- 
volved in  the  balance  consideration  for  the  indicated  layer. 

In  order  to  compute  F26  and  F610  at  each  gridpoint,  complete 

*     *     * 
listings  of  IR  net  fluxes  Fn  _  ,  F.  ,  F_  ,  associated  with  each  set  of 

10    6    2 

reference  cloud  amounts  (0,0),  (1,0),  (0,1),  (1,1)  have  been  computed 

* 

at  each  gridpoint  together  with  the  weighted  set  of  fluxes  F    (CL  , 

*  * 

CL  ) ,  F   (CL  ,  CL  )  ,  F   (CL  ,  CL  )  .   These  cloud-weighted  values  have 

£»  O       JL       4*       £*  JL       £* 

been  constructed  using  the  weighting  scheme  of  Eqs.  2-6,  3-1,  3-4  and 
3-5.   A  sample  gridpoint  printout  of  the  net  fluxes  and  flux  divergences 
has  been  included  in  Table  II. 

The  determination  of  the  heat  balance  by  layers  has  been  de- 
ferred to  Section  VI,  where  the  results  are  displayed  in  cross- 
sectional  form  for  the  2/3-CL  case. 

D.   STATISTICAL  RESULTS  AND  COMPARISONS 

1.   Net  Flux  F     (0,0) 

* 
A  statistical  test  of  the  F    (0,0)  numerical  results  was  con- 
ducted by  utilizing  a  linear  regression  program  for  relating  the 
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* 
predictand  F    (0,0)  against  the  simultaneous  values  of  X  =  B   and 

X  =  B  ve   as  predictors,  using  all  93  gridpoint  data  as  samples.   The 
regression  program  was  based  on  the  BMD02R  in  the  Biomedical  set  of 
programs  (Dixon  (1973)).   The  result  was  cast  in  the  form  of  the  well- 
known  Brunt  net  flux  equation  (1932): 

Fio*  =  B10(a+b'^)  (3"8) 


where 


4 
B    =  St*T  (3-9) 

10        10  v    ' 


is  the  Stefan-Boltzmann  blackbody  flux  at  the  surface  temperature  T 

and  e  is  the  surface  vapor  pressure  in  mb.   The  multiple  correlation 

*  i 
coefficient  R  (F    |X  ,X  )  was  higher  than  0.99  in  the  16  July  case, 

as  was  also  true  for  the  data  periods  of  Warner  (16  October  1973)  and 

Meyers  (16  April  1974) ,  (Table  III) . 

Since  F    may  be  defined  as 

Fio*  =  B10  "  Fd  (3"10) 

Eq.  3-10  can  then  be  solved  for  F_. 

d 

F^    =  B   r(l-a)  -  b^e"]  (3-11) 

d      10 

for  the  clear-sky  case.   It  must  be  recalled  that  the  computation  of 

F,  included  both  the  effects  of  water  vapor  and  of  CO^.   Since  there  is 
d  2 

no  direct  relationship  of  ve   upon  the  quasi-constant  absorber  mass  of 
CO  ,  the  radiative  effect  of  the  CO   in  Eq.  3-11  must  be  included  in 
the  term  involving  (1-a)  B   . 
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A  cross-seasonal  comparison  of  the  coefficients  1-a  and  -b  of 
Eq.  3-11  is  presented  in  Table  III,  comparing  the  clear-sky  results  of 
F   for  16  July  1975  with  those  of  Meyers  and  Warner  at  climatological  in- 
tervals of  three  months  earlier  and  later,  respectively. 

The  high  multiple  correlation  coefficients  indicate  that  an 

oceanic  version  of  the  Brunt  downward  flux  equation  is  valid  for  each 

tested  season.   F   depends  most  strongly  upon  B   while  only  slightly  on 

B   /e.   This  is  an  apparent  result  of  the  nearly  constant  relative 

humidity  over  the  large  oceanic  regimes.   Moreover,  the  weak  dependence 

of  F,  upon  B, „ve   is  most  evident  at  the  data-time  16  July  1974,  when 
d       10 

the  coefficient  -b  is  so  small  as  to  be  statistically  insignificant. 
This  result  is  due  primarily  to  the  relatively  weak  gradients  of  both 
/e  and  in  the  radiative  soundings  aloft  corresponding  to  16  July  1974, 

as  compared  with  the  adjacent  seasons. 

* 

2.   Modification  of  F     for  Cloudiness  CL  ,  CL 

A  second  statistical  test  was  performed  based  on  a  possible  rela- 

* 
tionship  between  the  cloudy-sky  and  the  clear-sky  cases  of  F    .   The 

*  * 

ratio  Y  =  F    (CL  ,CL  )/F    (0,0)  was  used  as  the  predictand  in  this 

case.   The  predictor  CL  was  the  total  opaque  cloud  cover  at  each  grid- 
point,  defined  as  (after  Quinn,  1971; and  Spaeth,  1975): 

CL  =  CL(1)  +  CL(2)  (l-CL(l)  )  (3-12) 

* 

This  regression  was  run  utilizing  F  (CL  ,  CL  )  in  the  numerator,  com- 
puted using  the  amounts  of  CL(1)  and  CL(2)  as  calculated  from  Smagorin- 
sky's  formulation  and  then  utilizing  F    (CL  ' ,  CL  * )  computed  from  the 
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2/3-CL  parameterization.   The  general  form  of  the  regression  formula 
required  is  shown  in  Eq.  3-13  below,  with  CL  given  in  Eq.  3-12 

F10*(CL)  =  F1()*(0,0)  (l-d(CL))  (3-13) 

A  cross-seasonal  comparison  of  the  coefficient  d   for  the  2/3-CL  and 
full-CL  cases  using  the  results  of  this  study  and  those  of  Meyers  and 
Warner,  and  of  the  resulting  multiple  correlation  coefficients  are  pre- 
sented in  Table  IV.   There  seems  to  be  little  clear-cut  variation  in 
the  statistics  of  Table  IV  between  the  results  of  adjacent  seasons, 
and  little  statistical  preference  between  the  full-CL  and  the  2/3-CL 
cases.   The  high  correlation  coefficients  of  Table  IV,  for  all  seasons 
considered,  indicate  the  general  capability  of  the  IR  radiative  model 
employed  in  this  study  to  account  for  the  effective  net  radiation  at 

the  surface. 

* 

The  mean  values  of  F    over  the  93  gridpoint  soundings  for 

16  July  1974,  as  applicable  in  Tables  III  and  IV  were  as  follows: 
Full-CL  Statistics  2/3-CL  Statistics 

F   *(CL)  =  .0774  ly  min*"1        F   *(CL')=  .0966  ly  min"1 


.  -1        ~   „,„  ^    .  -1 


F   *(0,0)=  .1449  ly  min  F   *(0,0)  =  .1449  ly  min 

CL  =  .5890  CL'=  .4306 

The  superior  bar  symbol  denotes  the  sample  mean.   These  results  indi- 
cate that  the  surface  net  flux  is  decreased  by  the  ratio  .4658  with 


CL  =  .5890,  and  by  the  ratio  .3333  with  CL'  =  .4306.   These  results 

indicate  that  the  downward  flux  F  (CL)  is  46.6  percent  larger  at  the 

d 

surface  than  F  (0,0)  when  full-CL  amounts  are  used  in  the  cloud  para- 
meterization model.   In  the  2/3-CL  parameterization  the  mean  downward 
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16  April  1974 


16  July  1974 


16  October  1973 


1-a 

.6557 

-b 

.0209 

F~ 
d 

.4527 

R 

.9969 

M 


.8129 

-.0094 

.5471 

.9950 


.6834 
.0161 
.4716 
.9973 


Table  III.  Cross-seasonal  comparisons  of  the  coefficients  1-a  and  -b 
from  Eq.  3-11,  as  well  as  mean  values  from  the  Brunt  equation  and  the 
resulting  multiple  correlation  coefficients  have  been  included. 


16  Apri 
2/3-CL 

.1  1974 
full-CL 

16  July 
2/3-CL 

1974 
full-CL 

16  October  1973 
2/3-CL    full-CL 

d 
R 

-.73 
.9916 

-.75 
.9945 

-.74 
.9942 

-.81 
.9920 

-.72 
.9974 

-.73 
.9972 

M 


Table  IV.   Cross-seasonal  comparison  of  the  coefficient  d   from 

Eq.  3-13  and  the  resulting  multiple  correlation  coefficients  for  both 

the  2/3-CL  and  the  full-CL  cases. 
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flux  at  the  surface  F.(CL')  is  increased  by  only  33.3  percent  when 

d 

compared  to  F  (0,0). 

E.   COMPARISONS  OF  FF2  WITH  SATELLITE  CLIMATOLOGY 

Comparison  was  made  of  computed-model  values  of  FF2  with  satellite 
measurements  of  total  long-wave  flux  to  space  (after  Raschke  et  al, 

1973)  for  the  NIMBUS  III  period  of  16-31  July  1969.   FF2  is  used  here 

* 

instead  of  F    since  the  downward  flux  at  k=2  is  not  measured  by  the 

satellite.   The  computation  of  FF2  was  previously  discussed.   FF2  for 
16  July  1974  was  computed  for  both  full-CL  and  2/3-CL  cases  at  all  93 
gridpoints.   These  FF2  values  were  then  interpolated  to  whole  multiples 
of  5°  latitude  along  each  of  the  four  meridians  and  the  results  were 
then  averaged  zonally  to  get  a  mean  latitudinal  distribution  valid  for 
each  5°  latitude  interval  considered  in  this  study. 

Table  V  shows  both  of  the  latitudinally- distributed  FF2  cloud- 
model  values  compared  with  those  extracted  from  Raschke  et  al  (1973) . 
Raschke' s  results  were  obtained  by  interpolation  from  charts  using  the 
same  oceanic  meridians  as  those  used  in  this  study.   Latitudinally  dis- 
tributed total  opaque  cloud  cover  CL  by  Eq.  3-12  are  also  listed  for 
both  full-CL  and  2/3-CL  cases  at  each  latitude.   In  the  bottom  line 
of  each  column  in  Table  V  is  listed  the  cosine-weighted  mean  of  each 
set  of  the  column  values  for  the  listed  latitude  range.   It  should  be 
noted  that  the  Raschke  results  do  not  correspond  to  specifically 
known  CL  values.   Table  V  suggests  that  the  model  results  for  FF2, 
both  full-CL  and  2/3-CL,  are  reasonably  close  to  the  corresponding 
values  derived  from  Raschke.   A  close  examination  of  this  table  shows 
that  in  the  regions  20S  to  5N  and  25N  to  4  5N  both  sets  of  model  results 
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RADIATIVE  STATISTICS  CLOUD  STATISTICS 


Lat 

FF2 
full-CL 

FF2 
2/3-CL 

F2(RAS) 

full-CL 
CL  Avg 

2/3-CL 

CL'  Avg 

20S 

.377 

.382 

.40 

.398 

.269 

15 

.386 

.390 

.42 

.456 

.306 

10 

.400 

.403 

.42 

.396 

.258 

5 

.378 

.384 

.39 

.616 

.498 

0 

.357 

.366 

.38 

.615 

.445 

5 

.341 

.354 

.36 

.737 

.534 

10 

.357 

.367 

.34 

.651 

.460 

15 

.369 

.377 

.35 

.622 

.430 

20 

.384 

.389 

.37 

.638 

.431 

25 

.390 

.393 

.40 

.559 

.375 

30 

.372 

.378 

.39 

.601 

.425 

35 

.378 

.382 

.38 

.555 

.383 

40 

.356 

.360 

.36 

.546 

.387 

45 

.333 

.338 

.36 

.671 

.566 

50 

.284 

.298 

.35 

.652 

.475 

55 

.210 

.229 

.33 

.948 

.819 

60 

.269 

.286 

.33 

.635 

.459 

65 

.196 

.219 

.33 

.862 

.590 

Wt 

Avg 

.352 

.360 

.374 

.552 

.405 

Table  V.   Comparison  of  zonally-averaged  longwave  flux  to  space,  FF2, 
as  found  by  this  study  for  both  full-CL  and  2/3-CL  cloud  cases  for 
16  July  1974,  and  by  Raschke  et  al  (1973)  based  upon  NIMBUS  III  measure- 
ments.  Also  included  are  composite  cloud  amount  fractions  used  in  the 
two  cloud  models  of  this  study.   Flux  values  in  ly  min"l. 
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are  slightly  less  than  given  by  Raschke,  while  in  the  region  ION  to 
20N  both  sets  of  values  are  slightly  higher  than  Raschke.   From  20S 
to  45N  the  agreement  between  FF2  and  F2(RAS)  is  quite  close.   From 
50N  and  higher  the  correspondence  to  F2(RAS)  is  not  as  good  and  perhaps 
this  is  due  to  the  high  amounts  of  cloud  cover  experienced  on  this 
day  in  this  region.   This  could  be  particularly  true  since  the  model 
calculations  are  for  a  single  day,  whereas   Raschke 's  calculations 
are  based  upon  16-day  averaged  values.   Additional  limitations  of  the 
comparisons  made  here  are  obvious,  when  it  is  recalled  that  between 
latitudes  20S-5S  and  between  60N-65N  there  are  fewer  than  four  meri- 
dional lines  available  for  computing  the  listed  zonal  values  in 
Table  V.   For  all  other  latitudinal-average  values,  four  meridional 
lines  were  used  in  the  averaging.   While  the  results  from  Table  V  are 
not  used  to  make  a  final  choice  of  cloud-cover  parameterization  for 
use  in  the  radiation  model,  there  is  evidence  of  slightly  closer  agree- 
ment of  FF2  with  F2(RAS)  using  the  2/3-CL  parameterization. 

* 

F.   COMPARISONS  OF  CROSS-SEASONAL  RESULTS  OF  F 

10 

* 
Table  VI  depicts  the  latitudinally  distributed  values  of  F    (CL') 

and  of  CL'  at  the  earth's  surface  obtained  using  the  2/3-CL  parameteri- 

* 
zation.   A  seasonal  comparison  of  the  model  computed  F    is  made  with 

data  obtained  from  Meyers  (1975)  for  the  16  April  1974  case  and 

from  Warner  (1974)  for  16  October  1973.   It  is  clearly  shown  that 

* 

F    (CL1)  is  a  decreasing  function  of  cloud  cover.   The  most  noteworthy 

* 
aspect  of  Table  VI  is  the  belt  of  high  values  of  F    for  16  July  1974 

in  the  latitudinal  band  20S  to  10S.   If  the  associated  CL  values  are 

examined,  it  is  noted  that  the  composite  cloud  values  were  very  low 
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16  April 

1974 

16  July 

1974 

16  October 

1973 

* 

* 

* 

Lat. 

F1Q  («•) 

CL' 

F1Q  (CL') 

CL' 

F10  (CL') 

CL' 

20S 

.0736 

.842 

.1509 

.269 

.0927 

.628 

15 

.0667 

.838 

.1328 

.306 

.0823 

.629 

10 

.0710 

.752 

.1346 

.258 

.0807 

.607 

5 

.0690 

.647 

.0935 

.498 

.0747 

.707 

0 

.0779 

.659 

.1014 

.445 

.0952 

.489 

5 

.0802 

.591 

.0945 

.534 

.1011 

.471 

10 

.0863 

.639 

.1022 

.460 

.1015 

.459 

15 

.0919 

.548 

.1033 

.430 

.0988 

.486 

20 

.0847 

.536 

.0941 

.431 

.1011 

.437 

25 

.0939 

.473 

.0983 

.375 

.1056 

.383 

30 

.0981 

.446 

.0918 

.425 

.1037 

.448 

35 

.1115 

.520 

.0924 

.383 

.1087 

.484 

40 

.1162 

.408 

.0854 

.387 

.1218 

.325 

45 

.0995 

.467 

.0705 

.566 

.1188 

.301 

50 

.0795 

.674 

.0637 

.475 

.1037 

.381 

55 

.1082 

.555 

.0377 

.819 

.0921 

.452 

60 

.0538 

.816 

.0852 

.459 

.0856 

.450 

65 

.0816 

.796 

.0613 

.590 

.0346 

.885 

Wt 

Avg. 

.0773 

.596 

.0848 

.405 

.0853 

.462 

Table  VI.   Cross-seasonal  comparison  of  the  zonally-averaged  net 


fluxes  at  the  surface  F    and  zonally-averaged  CL'  based  on  2/3- 
CL  parameterization  for  the  three  mid-seasonal  dates  16  April  197- 
16  July  1974,  and  16  October  1973.   Flux  values  are  in  ly  min~l. 
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* 
and  this  seems  to  be  consistent  with  the  high  values  for  F    .   Further- 
more, the  low  CL  values  would  appear  to  be  consistent  with  the  northward 

migration  of  the  ITCZ  during  mid-summer.   There  is  a  tendency  in  each 

* 
season  for  a  relative  maximum  value  of  F    to  be  located  in  the  sub- 
tropics  (latitudes  15-25IJ)  .  There  is  further  evidence  of  a  high  latitude 

* 
(55-65N)  minimum  F    associated  with  a  concentration  of  maximum  cloud 

cover  CL' .   The  analysis  does  seem  to  be  consistent  across  seasons.   Simi- 
lar conclusions  were  reached  for  the  full-CL  model,  although  they  were 

not  included  in  Table  VI.   In  general,  it  is  felt  that  the  results  of 

* 

the  F    analysis  yield  no  conclusive  results  in  regard  to  the  appropri- 
ateness of  the  choice  of  full-CL  or  the  2/3-CL  parameterization  models. 
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IV.   SOLAR  RADIATION 


A.   PARTITION  OF  SOLAR  INSOLATION 

The  solar  constant  assumed  in  this  study  at  level  k=0  (top  of  atmos- 
phere) was  2.00  ly  min    (Joseph,  1971) .   This  flux  was  depleted  by  4% 
to  account  for  the  attenuation  caused  by  oxygen  and  ozone  above  the  tropo- 
pause.   This  left  the  value  S  =  1.92  ly  min    at  level  k=2  to  be  used 
in  this  study  as  the  effective  solar  constant. 

Equation  (4-1)  was  then  utilized  to  compute  the  effective  solar  in- 
solation at  the  tropopause  (k=2)  as  follows: 

F(2)  =  S  [— ]~2  Cos  z  (4-1) 

r 
m 


where  S  =  effective  solar  constant  at  k=2 

Cos  z  =  cosine  of  the  zenith  angle  for  the  Julian  date  used 

r/r   =  ratio  of  the  actual  earth-sun  distance  to  the  mean  earth-sun 
distance  for  the  Julian  date  used  in  this  study. 

The  Smithsonian  Meteorological  Tables  (List,  1958)  gives  the  ratio  r/r 

m 

and  the  solar  declination  6  for  16  July  1974,  0000GMT,  as  these  values 
are  employed  in  Eq.  4-1.   These  values  are,  respectively: 


—  =  1.01644 
r 
m 

6  =  21. 5N  latitude 


6  is  used  in  evaluating  the  cosine  of  the  solar  zenith  angle,  given  by 
Cos  z  =  Sin  (J)  Sin  6  +  Cos  <J>  Cos  6  Cos  h  (4-2) 


46 


where   (j)   is  the  latitude  and  h   is  the  hour  angle  of  the  sun  relative 
to  the  meridional  data  lines.   For  example,  Fig.  1  makes  it  clear  that 
for 

line  1       h  =  55° 

line  2       h  =  10° 

lines  3,4    h  =  35° 
at  the  times  of  the  synoptic  charts  (0000GMT  and  1200GMT,  as  applicable) . 
Sin  (j)  was  computed  according  to  the  standard  polar  stereographic  projec- 
tion formula  applicable  to  the  FNWC  base  chart  given  by 

r  2  -  [(I-32)2  +  (J-32)2] 

Sin  <j>  =  — —  (4-3) 

r    +  [(1-32)   +  (J-32p] 

2 
where   r   =  973.752.   Thus  Sin  (J)  assumes  the  following  functional  form 
E 

in  terms  of  the  FNWC  grid-coordinate  I  (Fig.  1) : 

Lines  1,  3,  4 

0.   .    973.752  -2(32-I)2  ,  „  „ ,    .. 

Sm  (J)  =  s '—  (4-4  (a)) 

973.752  +2(32-1) 

Line  2 

2 

2 


_.   ,    973.752  -  (32-1) 

Sin  <p   = (4-4  (b)  ) 


973.752  +  (32-1) 
or,  conversely  I  is  given  in  terms  of   <j)   by 

Lines  1,  3,  4 

I  =  32  -  22.065  [-~|t^—]  (4-4  (c)) 

1  +  Sin  <p 

Line  2  I  =  32  -  31.205[1  ^^    J  (4-4  (d)  ) 

I  =  1,...,  25  for  Lines  1,  2 
I  =  9,...,  25  for  Line  3 
I  =  63 ,...,38  for  Line  4 
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For  lines  1,  2,  and  3  the  gridpoint  soundings  correspond  to  OOOOGMT  16 
July  when  solar  noon  occurs  at  the  180th  meridian.   Line  4  gridpoint 
soundings  correspond  to  1200GMT,   16  July,  when  solar  noon  was  at  the 
Greenwich  meridian. 

A  very  simple  partition  of  solar  insolation  was  utilized  in  this 
study  after  Joseph  (1971).   It  consisted  of  dividing  the  insolation  F(2) 
into  two  parts  at  level  k=2.   One  part  was  considered  to  include  all 
wavelengths   A  >  .9  urn  where  absorption  by  water  vapor  and  carbon 
dioxide  bands  are  the  most  prevalent  attenuation  processes  in  clear  air. 
This  part  of  the  solar  spectrum  was  termed  F(A)  energy  and  was  considered 
subject  to  water-vapor  absorption  but  not  to  Rayleigh   scattering.   For 
those  wavelengths  X  <  .9  ym,  absorption  of  the  solar  insolation  energy 
by  water  vapor  was  considered  negligible.   This  part  of  the  solar  insola- 
tion was  denoted  F(S)  suggestive  of  the  fact  that  it  was  subject  only  to 
Rayleigh  scattering  attenuation  in  clear  air.   The  two  partitions  are 
formulated  after  Joseph  (1971)  as  follows: 

F(A)  =  .349  F(2)  (4-5(a)) 

F(S)  =  .651  F(2)  (4-5  (b)) 

In  this  study,  the  introduction  of  two  cloud  decks  produced  cloud-reflec- 
tivity effects  upon  both  the  F (A)  and  F(S)  solar  energy  insolations. 
However,  in  the  clear  areas  around  any  gridpoint  the  absorption-attenua- 
tion only  applies  to  the  F(A)  insolation,  while  only  Rayleigh  scattering- 
attenuation   applies  to  the  F(S)  insolation. 

B.   DISPOSITION  OF  F(S)  INSOLATION 

In  the  disposition  of  the  F(S)  insolation,  Joseph  (1971)  determined 
that  Rayleigh  scattering  reflectance  to  space  by  clear  skies  (after 
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Coulson,  1959)  could  be  effectively  approximated  by  least  squares  in 
the  following  form 

a(R)  =  .085  +  .25074  [log  (~-   Sec  z)  ]  (4-6) 

o 

where   P   =  1013.25  mb.   In  Equation  4-6,  tt/P   =  1   in  view  of  the  fact 
o  ^  o 

that  the  mean  sea  level  pressure  7T   is  close  to  1000  mb.   Also 

Sec  z  =  (Cos  z) 

with  Cos  z  given  by  Equation  4-2. 

The  surface  albedo  a(G)   is  another  reflective  parameter  utilized 
in  this  study.   Over  oceanic  areas  the  following  formula  for  a(G) 
after  Gates  et  al  (1971) ,  was  utilized: 

a(G)  =  max  {.06,  .06  +  .54  (.7  -  Cos  z) }  .  (4-7) 

1.   Clear  Sky  Case 

In  the  clear  sky  (0,0)  case  the  F(S)  insolation  was  subjected 
to  Rayleigh  scattering-reflectance  a(R)   by  the  atmosphere,  as  well 
as  to  reflectance  0t(R)   by  the  earth's  surface.   Considering  the  likeli- 
hood of  a  succession  of  multiple  reflections  between  earth  and  atmosphere, 
the  F(S)  insolation  actually  penetrating  the  earth's  surface  after 
scattering  is  given  by 

IS10(0,0)  =  F(S)  [l-a(R)]  [l+a(R)a(G)  + (a(R)a(G))n 

+...]*(l-a(G))  (4-8(a)) 

that  is,  by 

IS10(0,0)  =  F(S)  [l-a(R)]  [l-a(G)]/[l-a(R)a(G)J  (4-8(b)) 
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2.   Cloudy-Sky  Cases 

In  the  three  cases  in  which  clouds  were  present,  F(S)  insola- 
tion absorbed  by  the  ground  at  each  gridpoint  was  (Arakawa,  1972) : 

IS10(1,1)  =  F(S)  (l-R(D)  (1-R(2))  (l-Ct(G)) 

*{l-[  (R(1)R(2)  +  R(2)a(G)  +  R(l)a(G) 

+  2R(l)R(2)a(G)]}_1  (4-9) 

As  indicated  by  the  notation  (1,1),  denoting  CL(1)  -   CL(2)  -   1.0, 
Eq.  4-9  is  the  formula  used  in  calculating  F(S)  insolation  absorbed  by 
the  ground  in  the  case  where  overcast  clouds  are  present  at  both  levels 
of  Fig.  2.   Also  in  Eq.  4-9,  constant  cloud-reflectance  values  were 
chosen,  namely  R(l)  =  .54  for  the  mid-level  clouds  between  k=4  and  6, 
and  R(2)  =  .66  for  the  low-level  clouds  between  k=8  and  9.   Both  cloud- 
reflectance  values  are  as  suggested  by  C.  D.  Rodgers  (1967) . 

For  all  the  other  cloud  cases,  the  following  changes  were  applied 
to  Eq.  4-9.   In  the  (1,0)  case  (CL(1)  =  1.0,  CL(2)  =  0.0),  the  desired 
result  is  obtained  by  setting  R(2)  =  0  in  (4-9) ,  from  which  it  follows 
that 

IS10(1,0)  =  F(S)  (l-R(l))  (l-a(G))/[l-R(l)a(G)  ]  (4-10) 

In  the  case  (0,1),  one  sets  R(l)  =  0.0  in  (4-9)  so  that  (4-9)  is  simpli- 
fied to 

IS10(0,1)  =  F(S)  (1-R(2))  (l-a(G))/[l-R(2)CX(G)]  (4-11) 

Note  that  with  a  cloud  overcast  present,  the  Rayleigh  clear-sky  scatter- 
ing cx(R)   does  not  appear  in  Eqs.  4-9,  4-10  or  4-11,  but  is  included 
in  the  cloud  reflectances  (Rogers,  1967) . 
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3.   Composite  F(S)  Insolation 

Equations  (4-8),  (4-9),  (4-10),  and  (4-11)  were  utilized  in  the 
computation  of  the  cloud-weighted  F(S)  insolation  penetrating  the 
earth's  surface  considering  the  areal-weights  of  the  cloud  combinations 
denoted  by  (0,0),  (1,1),  (1,0)  and  (0,1)  about  a  gridpoint.   The  resultant 
F(S)  insolation  penetrating  the  earth's  surface  denoted  by  IS10  is  there- 
fore expressible  as 

IS10(CL(1),  CL(2))  =  IS10(0,0)  W(0,0) 

+  IS10(1,1)  W(l,l) 

+  IS10(1,0)  W(1,0) 

+  IS10(0,1)  W(0,1)        (4-12) 

Here  the  weighting  factors  W(0,0),  W(l,l),  W(1,0)  and  W(0,1)  are  as  com- 
puted in  Eq.  2-6 (a,..., d)  respectively.   Note  that  Table  VII  lists  the 
radiative  computations  as  they  apply  to  the  disposition  of  F(S)  insola- 
tion.  The  four  individual  computations  of  IS10  as  applicable  to  the 
possible  overcast-clear  layer  cases  are  made  under  that  heading  while 
the  composite  value  of  IS10  appears  in  the  final  line  of  the  same  column 
of  Table  VII.   The  difference 

REFS  =  F(S)  -  IS10  (4-13) 

in  each  case  represents  F (S) -insolation  reflection  to  space  while 

STRAN  =  (4-14) 

(1  -  a(G))  K  ' 

has  been  computed  as  that  portion  of  the  F (S) -insolation  incident  at 
the  sea  surface  just  prior  to  absorption  by  the  surface.   Note  that 
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no  absorption  in  air  has  been  included  in  the  computations  of  Table  VII, 
and  that  the  only  absorption  permitted  is  that  implicit  on  IS10.   The 
final  line  under  each  column  of  Table  VII  depicts  the  composite  disposi- 
tions computed  by  means  of  the  weighting  scheme  of  Eq.  4-12  as  applicable 
to  the  case  2/3-CL. 


Table  VII.   A  sample  listing  of  F(S)  insolational  disposition  computed 
at  gridpoint  (2,2)  using  equations  listed  in  Sec.  IV. B.   Insolational 
values  in  ly  min~  . 


(CL  ,  CL  )     Weighting    IS10 


REFS 


STRAN 


(0,0) 

.8391 

.3203 

.1639 

.4117 

(1,0) 

.0000 

.1969 

.2873 

.2531 

(0,1) 

.1609 

.1501 

.3342 

.1929 

(1,1) 

.0000 

.110 

.3742 

.1414 

COMPOSITE  VALUES 

.2969 

.2217 

.3815 

C.   DISPOSITION  OF  F(A)  INSOLATION 

The  fractional  portion  of  the  solar  insolation  subject  to  absorp- 
tion by  atmospheric  water  vapor  and  carbon  dioxide  are  covered  in  the 
following  subsections. 

1.   Clear-Sky  Case  (0,0) 

The  Manabe-Moller  absorptivity  function  provided  the  necessary 
absorptivity  values  for  the  key  layers  in  this  case.  The  form  of  this 
absorptivity  function  is 


a(2,  k)  =  .271[U(2,  k)  Sec  z] 


.303 


(4-15) 


Here   a_  is  the  absorptivity  applied  to  the  pressure-scaled  water  vapor 
mass  between  levels  2  and  k  (Fig.  2)  along  the  zenith  slant-path  angle  z, 
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The  resultant  absorbed  insolational  energy  in  the  partular  layer  (2,6) 
is  then  given  by  the  Manabe-Moller  relation 

303 
A(2,6)  =  0„271F(A)  [U(2,6)  Sec  z]  *  (4-16) 

The  two  layers  of  interest  in  which  absorption  was  computed  were  (2,6) 
and  (2,10).  The  absorbed  insolation  in  the  layer  (6,10)  was  then  com- 
puted by 

A(6,10)  =  A(2,10)  -  A(2,6)  (4-17) 

Water  vapor  mass  above  level  2  was  assumed  negligible  in  the  F(A)  dis- 
position of  the  solar  insolation. 

By  subtracting  A (2, 10)  from  F (A) ,  the  direct  transmission  of 
F(A)  insolation  impinging  at  the  earth's  surface  was  determined.   The 
transmission  of  F(A)  insolation  is  then  further  reduced  by  the  transmis- 
sivity  of  (l-a(G)),  after  surface-reflectance,  which  leads  to  the  earth- 
absorbed  result 

303 
IA10(0,0)  =  F(A){l-.271[U(2,10)  Sec  z]  *    }(l-a(G))       (4-18) 

The  transmitted  energy  impinging  upon  the  earth  just  prior  to  absorption 
is 

TRANA(0,0)  =  IAl0(0,0)/[l-a(G)]  (4-19) 

2.   Cloudy  Cases 

In  order  to  compute  meaningful  dispositions  of  F(A)  insolation 
in  cloudy-sky  cases,  cloud  reflectivities  and  cloud  absorptivities 
after  C.  D.  Rodgers  (1967)  were  utilized.   The  reflectance-values  used 
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are  different  from  those  suggested  by  Rodgers  for  the  F(S)  wavelengths. 
The  cloud  reflectivities  used  here  are  RA(1)  =  .46  and  RA(2)  =  .50.   In 
this  case  there  are  also  cloud-absorptivities  to  be  considered.   These 
were  taken  as  A(l)  =  .20  and  A (2)  =  .30,  respectively.   In  the  following 
discussions  the  cloud  conditions  are  considered  totally  overcast  in  the 
layer  or  layers  indicated  when  the  notations  ((1,1),  (0,1),  and  (1,0)) 
are  utilized. 

Schematic  representation  of  the  computations  performed  in  the 
various  cases  ((1,1),  (0,1),  (1,0))  are  displayed  in  Fig.  3.   This  figure 
indicates  the  theoretical  parameters  required  in  computing  the  disposi- 
tion of  incoming  F (A) -insolation  from  level  k=2  to  the  earth's  surface 
(k  =   10) .   The  equations  which  relate  to  the  parameters  in  Fig.  3  and' 
the  other  cloud  configurations  are  listed  in  Appendix  B  of  Spaeth's 
study  (1975) . 

The  symbols  A24,  A46,  A68,  A89,  and  A910,  etc.,  represent  the 
insolation  absorbed  in  the  layers  concerned.   Symbols  F2,  F4,  F6,  F8  and 
F9,  etc.,  depict  the  streams  of  insolation  passing  through  the  indicated 
levelo   A  vertical  arrow  implies  the  direction  of  insolation  passage, 
i„e.,  i  denotes  downward  insolation,  t  upward-reflected  insolation,  and 
II  downward-reflected  insolation.   The  absorption  quantity  A(6,8)ll,  for 
example,  indicates  absorbed  energy  remaining  in  (6,8)  from  a  downward 
reflected  beam. 

Since  multiple  reflections  occur  between  the  earth's  surface 
and  a  cloud  base,  or  between  cloud-layers  in  this  model  (Fig.  3) ,  the 
insolational  amount  reaching  the  lowermost  reflecting  surface  is 
affected  by  the  downward  reflected  beams.   Additional  insolational 
amounts  remaining  after  more  than  two  reflections  are  negligible,  and 
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Figure  3.   Schematic  representation  of  F (A) -insolation  disposition 
in  the  case  of  two  overcast  layers  (cloud-case  (1,1)). 
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are  not  included  in  the  model  computations.   Note  that  insolation  re- 
flected upward  from  a  lower  interface,  either  cloud- top  or  ground,  to 
the  base  of  a  higher  cloud  deck  was  not  subjected  to  absorption  by  or 
transmission  through  the  cloud.   This  simplification  resulted  in  slight- 
ly reduced  F(A)  insolation-reflectance  to  space. 

A  final  point  to  note  is  that  there  are  three  contributions  to 
the  absorption  of  solar  insolation  within  a  layer  between  any  two  re- 
flecting surfaces.   In  Fig.  3,  these  contributions  are  identified  by 
arrows  which  indicate  the  portion  of  path  being  crossed. 

For  each  of  the  possible  overcast-cloud  combinations,  the 
transmitted  F(A)  insolation  arriving  at  the  earth's  surface  may  be  de- 
fined using  the  following  notation  (TRANA)  as 

TRANA(1,1)  =  FlOl  +  F10-H  (4-20) 

TRANA(1,0)  =  F104-  +  F10H  (4-21) 

TRANA  (0,1)  =  F104-  +  F1044  (4-22) 

where  the  right  side  parameters  of  (4-20) ,  (4-21) ,  and  (4-22)  are 
derived  in  Appendix  B  of  Spaeth  (1975) . 

The  F(A)  insolation  absorbed  by  the  earth  in  each  cloud-case 
considered  is  derived  from 

IA10(1,1)  =  F10I  (l-CX(G))  +  F1044  (4-23) 

IA10(1,0)  =  F104-  (l-a(G))  +  F104-I  (4-24) 

IA10(0,1)  =  F10I  (l-a(G))  +  F10II  (4-25) 

In  the  last  three  equations,  the  quantity  F1044  is  small  enough  in  each 
case,  that  no  further  reflections  from  the  earth  were  considered. 
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3.   Composite  F(A)  Layer-Absorptions  and  Surface- 
Absorption  Insolation 

As  has  been  previously  discussed,  the  standard  grid-area  weight- 
ing scheme  of  this  study  was  applied  to  obtain  composite  values  of  the 
absorbed  F(A)  insolation  in  key  layers  and  also  within  the  earth's 
surface.   The  weighting  factors  applied  to  the  corresponding  overcast- 
combination  absorption  quantities  provided  the  following  composite 
results: 

A26(CL(1)  ,CL(2))  -  A26(0,0)  W(0,0)  +  A26(l,0)  W(1,0) 
+  A26(0,l)  W(0,1) 
+  A26(l,l)  W(l,l)  (4-26) 

A610(CL(1) ,CL(2))  =A610(0,0)  W(0,0)  +A610(1,0)  W(1,0) 

+A610(0,1)  W(0,1) 
+A610(1,1)  W(l,l)  (4-27) 

IA10(CL(1)  ,CL(2)  )  =IA10(0,0)  W(0,0)  +  IA10(1,0)  W(1,0) 

+  IA10(0,1)  W(0,1) 
+IA10(1,1)  W(l,l)  (4-28) 

The  weighting  factors  W  (0,  0)  ,  . .  .  ,W  (1, 1)  were  first  listed  in  Eqs.  2-6 
(a,...,d).   Also  A26(0,0),  A610(0,0),  IA10(0,0)  are  given  in  each  clear- 
sky  case  (0,0)  about  each  gridpoint  by  Eqs.  4-16,  4-17  and  4-18, 
respectively. 

The  results  for  the  absorption  in  layers  (2,6),  (6,10)  and  at 
the  surface  are  shown  in  Table  VIII  for  the  four  cloud  cases  [(0,0), 
(1,0),  (0,1),  (1,1)]  and  for  the  composite  case  based  on  the  weighting 
factors  W(0,0) ,.. .,W(1,1) . 
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Table  VIII.   A  sample  listing  of  F(A)  disposition  parameters  computed 
at  gridpoint  (2,2)  using  equations  listed  in  Sec.  IV. C.   Values  listed 
are  in  ly  min~l. 


(CL  ,  CL  )    Weighting 


A26 


A610    AI10 


REFA    TRANA 


(0,0) 

.7587 

.0509 

.0670 

.1103 

.0314 

,1418 

(1,0) 

0.0 

.0821 

.0335 

.0424 

.1017 

.0538 

(0,1) 

.2413 

.0529 

.1215 

.0211 

.0641 

.0267 

(1,1) 

0.0 

.0821 

.0567 

.0102 

.1107 

.0129 

COMPOSITE 

VALUES 

.0514 

.0801 

.0888 

.0393 

.1260 

In  the  computational  scheme  indicated  by  the  entries  of  Table 
VIII,  the  reflected  F(A)  insolation  to  space  has  been  depicted  by  the 
symbol  REFA,  and  its  values  follow  from 


REFA  =  F(A)  -  A26  -  A610  -  IA10 


(4-29) 


whereas  the  TRANA  dispositions  are  given  by  Eqs .  4-19, ...  4-22,  respec- 
tively, or  by  its  weighted-mean  value  in  the  case  of  the  composites 
(last  line,  Table  VIII). 

4.   Absorptivity  (ABA)  by  Layers 

Here  the  (fractional)  absorptivity  as  well  as  the  actual  insola- 
tion values  absorbed  in  the  layers  are  considered.   In  the  computation 
of  absorptivity,  which  is  fractional  absorption,  the  total  undepleted 
insolation  at  the  top  (k=0)  is  used.   The  following  equation  was  utilized 
in  this  calculation: 


-2 
FADJ  =  2.00(r/r  )    Cos  z 
m 


(4-30) 


The  absorptivity  of  the  troposphere  ABA  was  computed  from  the  ratio  of 
the  insolation  absorbed  in  the  troposphere  to  the  insolation  incident 
at  the  top  of  the  atmosphere  rather  than  at  k=2: 


*»  .  A<2.6)^M6,10)  (4-31, 


D.  ALBEDO  (ALB)  OF  THE  EARTH-TROPOSPHERE  SYSTEM 

In  considering  the  planetary  albedo,  the  reflected  insolation  of 
the  earth-troposphere  system  in  both  the  F(A)  and  F(S)  insolational 
regions  must  be  recalled  by  the  program.   Thus  REF  is  computed  at  each 
gridpoint  as  the  sum  of  the  reflected  insolation  energy  in  F(A)  (denoted 
REFA  in  Eq.  4-29)  and  the  reflected  part  of  F(S),  previously  denoted 
REFS  in  Eq.  4-13: 

REF  =  REFS  +  REFA  (4-32) 

Finally  the  planetary  albedo  is  related  to  FADJ  through 

"*  "  5S5J  (4"33) 

E.  COMPOSITE  ABSORPTIVITY  (ABG)  BY  THE  EARTH-SURFACE;  COMPOSITE 
ATMOSPHERIC  TRANSMISSIVITY  (ATRAN) 

1.   Absorptivity  (ABG)  of  Earth 

By  summing  the  weighted  values  of  F(S)  and  F(A)  portions  of  the 

incoming  insolation  entering  the  earth,  the  total  insolation  absorbed 

at  the  earth's  surface  was  computed.   This  quantity  when  divided  by  the 

extraterrestrial  insolation  gave  the  fractional  absorptivity  (ABG)  of  the 
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earth's  surface.   The  equation  for  ABG  was 

IA10   +   IS10  ..    ,,, 

880  - ¥K5T~  ",-34) 

where  IA10,  IS10,  and  FADJ  were  defined  previously  by  Eqs.  4-28,  4-12, 
and  4-30,  respectively. 

2.   Transmissivity  (ATRAN)  of  the  Troposphere 

Also  computed  with  the  total  insolational  energy  TRAN,  incident 
at  the  earth's  surface  just  before  absorption  by  the  surface.   This  cal- 
culation is  given  by 

TRAN  =  TRANA  +  STRAN  (4-35) 

Here   STRAN  =  (IS10/ (1-a (G) )   as  previously  defined  in  Eq.  4-14,  and 
a(G)   is  given  in  Eq.  4-7.   TRANA  has  also  been  defined  as  the  weighted 
value  of  TRANA(0,0),  TRANA (1,1),  TRANA (1,0),  and  TRANA(0,1)  given  by 
Eqs.  4-19,. o., 22.   Note  also  in  justification  of  STRAN  that  the  four 
cases  for  IS10  of  Eqs.  4-8,..., 11  each  have  the  common  factor  (l-a(G)) 
in  the  numerator  and  therefore  each  transmitted  F (S) -insolational  com- 
ponent available  at  the  earth  just  before  absorption  needs  only  to  be 
divided  by  (l-a(G)).   TRAN  may  thus  be  viewed  as  the  total  insolational 
energy  incident  at  a  pyrheliometer  located  at  earth.   The  (fractional) 
transmissivity  of  the  troposphere  (ATRAN)  is  then  computed  from 

ATRAN  =  TRAN/FADJ  (4-36) 

Note  finally  that  the  major  dispositions  of  the  total  insolation  at 
the  indicated  map  times  have  now  been  identified  by  the  fractional 
values,  ALB,  ABA  or  ABG,  and  ATRAN,  each  of  which  is  a  fractional  value 
representing  the  reflectivity  (albedo) ,  absorptivity  of  air  or  earth, 
or  atmospheric  transmissivity  as  the  case  may  be. 
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3 .   Computational  Check 

The  computational  scheme  utilized  in  this  model  was  checked  con- 
tinuously by  summing  the  fractional  values  ALB,  ABA,  and  ABG  at  each 
gridpoint.   The  value  in  each  case  must  total  .96,  since  as  previously 
noted  the  attenuation  of  solar  insolation  was  taken  as  four  percent  as 
it  passed  through  the  stratosphere. 

F.   STATISTICAL  ANALYSIS 

In  order  to  substantiate  some  of  the  computations  performed  in  this 
section,  several  of  the  most  important  items  computed  were  statistically 
analyzed  using  linear  regression  computer  programs  from  the  BMD  set  of 
statistical  programs  (Dixon,  1973). 

lo   Clear  Sky  Cases 

Using  ALB  (0,0,  z)  which  is  the  clear  sky  case  of  albedo,  as  the 
predictand,  and  log   Sec  z  as  the  predictor,  the  following  best-fit 
equation  resulted 


ALB(0,0,z)  =  .07856  +  .38689  log   Sec  z 

R   =  .9618  (4-37) 

M 


where   R   represents  the  multiple  correlation  coefficient. 

This  result  was  expected  over  the  ocean  where  a(G)   and  a (R) 
involved  logarithmic  dependence  on  Sec  z.   The  sample-average  values 
of  ALB  and  log   Sec  z  were 


ALB(0,0)  =  .1295  log   Sec  z  =  .1316 
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Other  clear-sky  regression  tests  made  use  of  the  water-vapor  mass  path 
length  (M)  defined  by 

1/2 
M  =  (U  Sec  z)     log   (U  Sec  z)  (4-38) 

This  parameterization  of  water  vapor  mass  is  similar  to  that  developed 
by  Hanson  (1971) ,  who  used  a  similar  M  as  a  predictor  in  his  empirical 
formulations  of  ABA  for  both  clear  and  partly  cloudy  sky  cases.   In  this 
study  the  two  regressions  attempted  using  M  as  a  predictor  were  with 
ABA (0,0, M)  and  ATRAN(0,0,M)   as  predictands.   The  best-fit  equations 
resulting  were 

ABA(0,0,M)  =  .1028  +  .0281  M 

R  =  .9922  (4-39) 

M 

ATRAN(0,0,M)  =  .7976  -  .0338  M 

R,  =  .8796  (4-40) 

M 

The  absorptivity  in  the  form  given  by  (4-39)  agrees  closely  in  form 
with  that  obtained  by  Spaeth  based  on  16  January  1974  data  and  with 
that  of  Hanson  (1971)  based  on  pyrheliometric  data.   The  means  of 
ABA (0,0, M)  and  ATRAN(0,0,M)  for  the  original  93  gridpoint  sample  were 


ABA(0,0,M)  =  .1433 


ATRAN(0,0,M)  =  .7488 


—  -2  1/2 

M  =1.44376  (gm  cm   )  / 


2.   Statistical  Relationships  Between  ALB,  ABA  and  ATRAN 
in  the  Cloudy  and  Clear  Sky  Cases 

In  addition  to  computing  the  clear-sky  value  of  ALB,  it  was  also 

possible  to  compute  the  composite  cloud  value  of  ALB.   A  regression  was 
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formed  showing  the  relationship  between  the  ratio  ALB(CL(1),  CL(2))/ 
ALB (0,0)  as  the  predictand  and  the  total  opaque  cloud  cover (CL)  as  the 
predictor.   CL  in  this  model  is  specified  by  Eq.  3-12,  and  gives  a  good 
approximation  of  the  effective  total  cloud  cover  by  the  two  layers  of 
cloud  amounts  CL  and  CL  by  Eq.  2-4  used  in  this  model.   As  explained 
earlier  in  this  study,  a  one-third  reduction  of  both  CL(1)  and  CL(2) 
after  initial  determination  by  Eq.  2-4  was  also  tried  in  all  computa- 
tions involving  cloud  amounts.   This  test  represented  an  attempt   at 
tuning  the  cloud  model  for  radiative  calculations,  and  gave  results  in 
closer  agreement  with  the  latest  literature  on  satellite  reflectances. 
Thus  all  subsequent  regressions  of  ALB (CL) /ALB (0,0)  upon  CL,  computed 
after  Eq.  3-12,  will  include  best-fit  equations  formulated  using  a  CL 
based  upon  the  2/3-CL  model  (Eq.  2-5)  and  the  full-CL  model  (Eq.  2-4) . 

A  small  letter  subscript  'a'  in  the  equation  number  or  symbol 
will  indicate  those  computations  using  2/3-CL  amounts  while  a  small 
letter  subscript  'b'  will  indicate  those  utilizing  the  full-CL  model. 

The  first  regression  equations  tested  were  those  for  ALB  (CL') 

a 

and  ALB  (CL) for  the  respective  cases  of  2/3-CL  and  full-CL  cases.   The 
b 

results  are  given  in  (4-41a,b) .   The  symbol  R  once  again  signifies 
the  multiple  correlation  of  the  statistical  regressions. 

ALB  (CL(1)CL(2))  =  ALB(0,0)[1  +  2.9915  CL  +  .5873  CL2] 
a 

R  =  .9858  (4-41(a)) 

M 

ALBl_(CL(l)CL(2)  )    =  ALB(0,0)  [1   +    2.7654   CL  +    .8308   CL2] 
b 

R      =    .9846  (4-41 (b) ) 

M 
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The  sample  means  of  the  values  in  Eqs.  4-41 (a, b)  are 
Eq.  4-41a  Eq.  4-41b 


ALB(CL(1) ,CL(2))  =  .3163       ALB (CL (1) ,CL (2) )  -  .3893 
a  b 


ALB(OfO)  =  .1295      ALB(0,0)  =  .1295 

CL  =  .4320      CL  =  .5890 

a  b 

The  high  value  of  the  correlation  coefficient  in  Eq.  4-41  shows  strong 
dependence  of  ALB  upon  CL  in  either  cloud  parameterization  case.    Also 
note  that  ALB (CL (1) ,CL(2) )  values  computed  here  are  clearly  higher  than 
reported  by  Raschke  et  al  (1973),  where  ALB  =  .26  for  essentially  the 
same  set  of  gridpoints  during  the  NIMBUS  III,  period  16-31  July  1969. 
The  2/3-CL  cloud  parameterization  gives  albedo  results  which  are  closer 
to  Raschke 's  results,  suggesting  that  the  full-CL  amounts  of  Smagorinsky 
(Eq.  2-4)  are  too  high  for  radiational  computations.   A  second  regres- 
sion was  developed  between  ABA (CL) /ABA (0, 0)  versus  CL.   The  following 
best-fit  equation  resulted: 

2 
ABA  (CL(1) ,CL(2))  =  ABA (0, 0) [1+ .4667CL  -  .0148CL  ] 
a 

R   =  .9565  (4-42(a)) 

2 
ABA  (CL(1) ,CL(2) )  =  ABA (0, 0)  [1+. 5606CL  -  .1238CL  ] 

R  =    .9597  (4-42 (b) ) 

M 

Mean  statistics  for  these  regression  cases  are 
Eq.  4-42a  Eq.  4-42b 


ABA  (CL(1),CL(2))  =  .1721      ABA  (CL (1) ,CL (2)  =  .1827 


ABA(0,0)  =  .1433      ABA(0,0)  =  .1433 

CL  =  .4320     CL,  =  .5890 
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From  the  means  of  (4-42) ,  it  is  seen  that  the  model  specifies  an  in- 
crease in  atmospheric  (solar)  absorptivity  with  increasing  cloud  cover. 
This  result  is  in  agreement  with  studies  of  Warner  (1974)  and  Spaeth 
(1975)  and  with  the  earlier  study  of  London  (1957) . 

An  analogous  statistical  regression  was  then  developed  for  the 
cloudy-sky  transmissivity  relative  to  the  clear-sky  transmissitivity . 
The  resulting  regressions  were 

ATRAN  (CL(1) ,CL(2) )  =  ATRAN (0, 0) [1- .6522CL] 
a 

R  =  .9982  '  (4-43 (a)) 

M 

ATRAN  (CL(1) ,CL(2) )  =  ATRAN (0, 0) [1- .6647CL] 
t> 

RM  =  .9977  (4-43 (b) ) 

M 

and  the  corresponding  mean  values  of  the  transmissivities  are  listed 
below: 

Eq.  4-43a  means  Eq.  4-43b  means 


ATRAN  (CL(1) ,CL(2) )  =  .5389    ATRAN  (CL (1) ,CL (2) )  =  .4576 
3.  Id 


ATRAN(0,0)  =  .7488    ATRAN(0,0)  =  .7488 

The  parameter  ATRAN (CL' ) /ATRAN (0, 0)  given  by  (4-43a)  is  in  very  good 
agreement  with  the  statistical  result  of  Savino-Angstrom,  which  was 

Q  =  Q  (1  -  .66CL) 

for  the  tropical  ocean  areas.   The  coefficient  of  CL  in  the  last  result 
is  known  to  be  a  decreasing  function  of  latitude  (Budyko,  1958) . 

Similarly,  the  parameters  ALB (CL1 ) /ALB (0, 0)  and  ABA (CL' ) /ABA (0, 0) 
behave  consistently  with  respect  to  increasing  CL,  with  somewhat  more 


realistic  mean  values  for  ALB(CL')  and  ABA(CL')  than  for  the  corresponding 
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full-CL  parameterization.   Hence,  in  the  mean,  the  statistical  analyses 
performed  seem  to  show  agreement  with  observational  results  of  other 
investigators.   Comparative  computations  of  the  atmospheric  transmissi- 
vity,  particularly  in  the  2/3-CL  case  appear  to  give  better  agreement 
with  ground-based  pyrheliometry . 

A  major  test  still  to  be  made  is  the  comparison  of  the  lati- 
tudinal-distribution of  ALB  (CL)  and  ALB,  (CL)  with  that  of  the  satellite 

a  b 

albedo  of  Raschke  et  al  (1973) . 

G.   ALBEDO  COMPARISONS  WITH  PUBLISHED  RESULTS 

The  tropospheric  albedo  computations  of  the  solar-insolation  model 
of  this  section  have  been  presented  for  both  full-CL  and  2/3-CL  cases, 
respectively.   These  albedos  have  been  interpolated  to  whole  multiples 
of  5°  latitude  between  20S  and  65N.   These  computations  were  made  for 
each  of  the  four  oceanic  meridians  and  the  resulting  albedos  averaged 
across  the  meridians  are  presented  as  a  function  of  latitude  in  Table 
IX  for  both  the  full-CL  and  2/3-CL  cases,  where  comparison  is  also  made 
with  satellite  albedos  of  Raschke  et  al  (1973) .   For  sake  of  consistency 
the  Raschke  albedo  data  have  been  extracted  from  atlas  presentations 
for  the  period  16-31  July  1969  to  5°  latitude  grids,  on  the  same  four 
meridians,  before  applying  the  zonal-averaging  process.   The  total 
opaque  cloud-cover  amounts  for  the  two  cases  are  listed  under  the 
headings  "full-CL"  and  "2/3-CL"  corresponding  to  the  composite  cloud- 
cover  formula,  Eq.  3-12. 

The  primary  results  to  be  noted  in  Table  IX  are  the  quasi-constant 
values  of  ALB(RAS)  between  20S,...,35N,  with  individual  latitudinal 
values  not  greatly  different  from  the  tropical-subtropical  mean  of 
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Table  IX.   Comparison  of  zonally-averaged  albedos  as  computed  by  this 
study  for  both  full-CL  and  2/3-CL  models  for  16  July  1974  with  those 
of  Raschke  et  al  (1973) ,  and  the  corresponding  values  of  full-CL  and 
of  2/3-CL,  each  as  a  function  of  latitude.   The  bottom  line  lists  the 
cosine-weighted  average  of  the  parameter  listed  in  the  tabular  column 
above. 

Albedo    Albedo    Albedo     Full-CL    2/3-CL 
Lat.     full-CL     2/3-CL     Raschke Avgc Avg, 


2  OS 

.335 

.296 

.20 

.398 

.269 

15 

.347 

.296 

.18 

.456 

.306 

10 

.320 

.268 

.16 

.396 

.258 

5 

.392 

.345 

.26 

.616 

.498 

0 

.372 

.303 

.24 

.615 

.445 

5 

.404 

.324 

.21 

.737 

.534 

10 

.368 

.292 

.25 

.651 

.460 

15 

.359 

.282 

.24 

.622 

.430 

20 

.362 

.281 

.22 

.638 

.431 

25 

.330 

.259 

.21 

.559 

c375 

30 

.337 

.272 

.21 

.601 

.425 

35 

.325 

.259 

.24 

.555 

.383 

40 

.317 

.257 

.25 

.546 

.387 

45 

.365 

.341 

.32 

.671 

.566 

50 

.391 

.309 

.34 

.652 

„475 

55 

.536 

.463 

.37 

.948 

.819 

60 

.356 

.288 

.40 

.635 

.459 

65 

.413 

.323 

.39 

.862 

.590 

Wt.Avg.   .364        .302        .25        .552       .405 
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ALB(RAS)  =  .22.   For  most  of  this  latitude  range  the  computed  values  of 
ALB(2/3-CL)  exceed  those  observed  by  Raschke  (Table  IX)  by  the  frac- 
tional amounts  .05  to  .10.   From  40N  poleward  ALB(RAS)  and  ALB(2/3-CL) 
are  in  fairly  close  agreement,  with  individual  values  close  to  the  local 
average  of  0.35.   Note  that  at  latitudes  60,  65N,  the  ALB(RAS)  results 
are  slightly  greater  than  ALB (2/3-CL)  which  may  be  due  to  an  effect 
associated  with  the  large  mean  zenith  angle  at  these  latitudes. 

As  noted  earlier  in  Table  V,  at  latitude  55N  there  is  a  localized 
maximum  of  CL'  associated  with  the  largest  computed  value  of  ALB (2/3-CL) 
=  .463.   The  later  value  is  not  verified  by  the  two-week  average  obser- 
vation of  ALB(RAS)  =  .37.   Overall,  the  computed  values  of  ALB(CL) 
exceed  those  of  ALB(RAS),  latitude-by-latitude,  in  much  the  same  degree 
as  indicated  by  cosine-weighted  means  of  Table  IX,  which  are  also 
listed  below: 


ALB(full-CL)  -  .36 


ALB (2/3-CL)   =  .30 


ALB(RAS)      =  .25 

To  summarize,  there  is  a  well-defined  preference  for  the  2/3-CL 
over  the  full-CL  model  in  giving  albedos  closer  to  those  observed  by 
satellite  for  the  same  general  period  (16-31  July  1969)  as  reported  by 
Raschke  et  al  (1973).   An  area  of  difficulty  even  with  the  2/3-CL 
model  occurs  in  tropical  and  subtropical  latitudes  where  Von  der  Haar 
and  Hanson  (1969)  have  reported  the  reality  of  smaller  than  suspected 
cloud-covers  and  of  smaller  values  of  the  resultant  global-albedo.   At 
these  latitudes  (tropical  and  subtropical)  it  seems  doubtful  that  the 
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large-scale  CL' -parameterization,  following  Smagorinsky,  is  as  applicable 
as  in  higher  latitudes.   Rather  some  method  of  specifying  the  distribu- 
tion (at  gridpoints)  of  isolated  cumulus  towers  may  be  appropriate. 
Finally,  it  is  questionable  whether  vertically-structured  convective 
cloud  elements  will  have  the  high  reflective  capability  attributed  here 
to  clouds  which  are  depicted  as  existing  in  horizontal  layers  with  a  con- 
stant reflectivity  .5  <  R  <  .65  as  specified  in  the  radiation  model.   It 
may  well  be  that  the  vertical  shafts  between  convective  elements  may  be 
a  more  efficient  focus  of  solar  energy  to  the  earth  than  is  perceived  by 
ascribing  a  constant  cloud  albedo  to  the  cloud-layer.   An  alternative 
approach  would  be  to  reduce  the  two-cloud  layering  in  the  tropics  to  a 
single  layer  of  clouds. 

At  any  rate  with  the  two-layered  cloud-model  considered  in  this 
study,  the  only  feasible  selection  from  the  results  deduced  here  is  the 
2/3-CL  parameterization.   This  selection  is  carried  over  into  the  compu- 
tations of  Section  VI.   A  further  test  of  the  cloud-parameterization 
effect  on  the  computation  of  net  radiation  of  the  earth-atmosphere 
system,  by  comparison  with  that  observed  by  satellite,  is  made  in  Sec- 
tion VII.  D. 
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V.   SENSIBLE  AND  LATENT  HEAT  TRANSPORT 
AT  THE  SEA-AIR  INTERFACE 

A.  GENERAL  PURPOSE 

The  model  used  here  to  describe  the  turbulent  transports  of  sensible 
and  latent  heat  across  the  sea-air  interface  was  basically  that  already 
in  use  in  the  operational  FNWC  primitive  equation  model.   This  model  has 
been  discussed  in  detail  by  Kesel  and  Winninghoff  (1972)  and  by  Kaitala 
(1974) .   The  primary  purpose  of  the  adaptation  of  the  turbulent  flux 
model  in  this  study  was  to  test  the  magnitudes  and  directions  of  these 
transfer  rates  in  comparison  with  those  of  the  radiational  transfer  model 
which  has  been  discussed  here  in  Sections  III  and  IV.   The  final  purpose 
of  the  inclusion  of  the  turbulent  flux  model  was  to  conduct  a  diagnostic 
computation  of  the  heat  balance  within  the  oceanic  and  atmospheric 
layers  for  16  July  1974  under  the  combined  operation  of  the  radiational 
and  turbulent  flux  models. 

B.  WINDSPEED  COMPUTATION  IN  THE  TURBULENT  FLUX  MODEL 

At  all  of  the  93  gridpoints  tested,  the  geostrophic  windspeed  V  was 
computed  by  use  of  the  standard  formula 

g    f  2(d/m) 
with  f  =  (ii£E_j  Sin  <j>  (5-2) 

1.86603 
m 


1+Sin  <f) 


and  AZ  =   /  (Z  -  Z2)2  +  (Z  -  Z±)z  (5-3) 
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In  Eq.  5-1,  d  =  381.0  km  is  the  nominal  FNWC  grid-spacing  (true  at  60N) 

and  ( — )  is  the  true-earth  spacing  used  in  the  centered-dif f erence  com- 
m 

putations  of  the  contour  gradient  (cf.,  Fig.  4). 

In  (5-1)  and  (5-3)   AZ   is  the  contour-difference,  centered  in  most 
cases  on  the  gridpoint  under  consideration.   However,  at  the  map  edge 
gridpoints  of  the  meridians,  it  was  only  possible  to  compute  a  forward- 
differenced  version  of  V   (cf.,  Fig.  4).   Thus  at  edge-point  "0",  Eq. 

d 
(5-4)  below  was  used.   The  spacing  —  in  Eq.  5-4  is  then  taken  as  one- 
half  that  in  (5-1) ,  so  that 

g  AZ  ._  .. 

g   f  (d/m) 

with 


AZ  =  /(Z4-  ZQ)2  +  (Z3-  ZQ)2  (5-5) 


All  geostrophic  windspeeds  were  then  converted  to  surface  windspeeds  us- 
ing the  following  empiricism  due  to  Langlois  and  Kwok  (1969) 

V  =  .8  V  +2.2  (5-6) 

s        g 

where  V  =  surface  windspeed  (mps)  and  V   is  the  1000  mb  geostrophic 
windspeed.   The  factor  2Q2  in  (5-6)  is  an  empiricism  introduced  to  in- 
clude the  effects  of  gustiness  upon  the  measured  (steady-state)  wind- 
speed. 

In  the  interval  of  +_  10  degrees  of  latitude,  the  value  of  the 

-4    -1 
Coriolis  parameter  was  arbitrarily  set  at  f  =  +_  .25  x  10   sec   ,  follow- 
ing Kaitala  (1974),  while  for  |<f)|  >  10°  ,  the  actual  value  of  f  for  the 
particular  latitude  was  used. 
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k  . 


•&*" 


Figure  4.   Section  of  FNWC  polar  stereographic  meridian  illustrating 
the  method  of  obtaining  contour  gradients  in  the  vicinity  of  a  grid- 
point.   The  values  of  contour-heights  correspond  to  the  1000  mb  D- 
values  in  the  neighborhood  of  the  gridpoint. 
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C.   SENSIBLE  HEAT  TRANSPORT 

The  details  of  the  adaptation  of  the  FNWC  model  for  sensible  heat 
transport  H„  utilized  in  this  study  follow  the  description  of  Spaeth 
(1975) .   The  formulation  of  H„  is  given  below  in  terms  of  bulk  para- 
meters of  the  surface  and  of  the  planetary  boundary  layer  (10,8)  of 
Fig.  2; 


Hr  "  pio  CP  cd  VTg  -  V 

"io  ■  pio(Rd  V1  (5"7) 


where  C      =    .239   cal   gm      (deg.    K) 

3 
and  C  =  1.4  x  10   (after  Weiler  and  Burling,  1967) 

The  equation  utilized  in  this  study  to  compute  T  may  be  shown  to  be 

(Spaeth,  1975) 

2KB 

C   V   T   +  —. ■ 2  K  Y 

T  =   D   S   9 ^ S.  (5-8) 

S+c   v 
AZ     D   s 


where 


e    =T    <I2oo  -286  (5-9) 

9     9   900  K         ' 


* 
k    = K    e  _  9  (5-10) 

9   10 

Az  =  V  zio  =  f    *  (ln  T]  (5-n) 


73 


*       * 

In  5-7 (c) ,  the  turbulent  parameters  K  and  a  have  been  determined 

empirically  by  FNWC 

*  •  ,«5   2    "I 

K  =  10  cm  sec 

a*  =  5  x  104cm(deg.K)~1 


T  -  T  in  Eq.  5-7  is  the  temperature  lapse  for  the  thin  layer  between 
g    x 

the  surface  and  the  top  of  the  constant  flux  layer  (cf.,  Fig.  5).   The 

formula  for  T  deduced  here  required  a  constant  value  of  sensible  heat 
x 

Hp  within  the  constant  flux  layer  (CFL)  at  each  gridpoint  with  Hr  then 
decreasing  linearly  with  pressure  to  H„  =  0  at  k=8,  and  then  remaining 
zero  above  k=8.   This  linear  decrease  of  Hp  with  pressure  is  compatible 
with  the  condition  of  constant  convergence  of  the  sensible  heat  flux  per 
gram  in  the  turbulent  boundary  layer  (x,8)0   Thus,  the  constant  value 
of  the  convergence  of  sensible  heat  as  given  by  Eq.  5-7  within  the  layer 
(8,10)  may  then  be  considered  to  be  applicable  at  level  k=9  as  depicted 
in  Fig.  5(a),  after  Spaeth  (1975). 

After  computing  T  by  Eq.  5-8,  Hp  was  then  computed  at  each  grid- 
point  using  the  bulk  transfer  equation  5-7.   The  gridpoint  values  of  H„ 
are  used  later  in  calculations  of  heat  balance  in  Sections  VI  and  VIII. 

D.   EVAPORATIVE  HEAT  TRANSPORT 

Similar  to  the  case  of  sensible  heat  convergence,  the  latent  heat 
transport  by  turbulence  was  assumed  to  be  subject  to  constant  evapora- 
tive flux  convergence  throughout  the  layer  800-1000  mb.   The  total 
amount  of  latent  heat  removed  by  evaporation  from  the  ocean  surface  was 
essentially  as  modeled  by  Kaitala  (1974)  using  bulk  transfer  theory. 
It  should  be  recognized,  however,  that  the  condensation  and  consequent 
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realized  latent  heat  may  occur  at  arbitrary  levels,  i.e.,  in  general  at 
levels  higher  than  k=8.   Thus,  in  following  the  FNWC  model,  we  have 
not  introduced  E  as  a  heating  rate  at  level  k=9. 

The  basis  of  the  turbulent  latent  heat  model  adopted  is  shown  in 
Fig.  5(b) ,  where 

E(900)  =  1/2  E  (Bulk)  (5~12) 

and  where  the  surface  layer  bulk  transfer  of  latent  heat  is  given  by 


E  =   E(Bulk)  =  L     C   V   (q  -q  )  (5-13) 

p    D   s   ^s  ^x 


Evaporation  at  level  k=9  may  also  be  computed  using 


-  <v  v 
E<9oo,-PloK(v^r 


where  K   is  given  in  Eq.  5-10.   In  (5-13)  and  (5-14) ,  L  is  the  latent 
heat  of  vaporization. 


L  =  596.73  -  .601  T 

in  cal(gm)   .   In  Eqs.  5-13,  5-14  and  5-15  q   and  q   are  the  mixing 

ratios  at  the  levels  k=9  and  at  x,  respectively,  whereas  q   is  the 

s 

saturation  mixing  ratio  at  the  surface  temperature  T  .   The  factor  1/2 
which  appears  in  the  right  side  of  (5-12)  reemphasizes  the  concept  of 
constant  latent-heat  convergence  in  the  layer  (x,8),  (cf.,  Fig.  5(b)). 
The  solution  for  q   is  obtained  by  combining  Eqs.  (5-12),  (5-13), 


and  (5-14)  in  the  form 


cd  vs  q10  +  2  h  % 

q  =  5 (5-15) 

X      C   V   +  2  ■£- 
D   s      AZ 


where  AZ  =  Z   -  Z    „ 
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Application  of  Eqs.  5-13  and  5-15  makes  possible  the  computation  of 

E  at  each  qridpoint.   As  in  the  case  for  T  ,  the  factor  of  2K  of  Eq. 

x 

5-15  was  taken  to  conform  to  the  value  K  as  used  by  Kaitala  (1974)  so 
that  the  turbulent  transport  results  obtained  here  will  be  compatible 
with  those  computed  operationally  by  the  FNWC  primitive  equation  model. 

As  already  noted,  this  was  accomplished  in  this  study  by  usinq  a  value 

5   2-1 
of  2K*  =  2  x  10  cm  sec   as  equivalent  to  the  K*-value  of  Kaitala  (1974). 

Thus  computation  of  E  by  Eqs.  5-13  and  5-15  was  equivalent  to  that  usinq 
Kaitala' s  scheme  (1974) ,  the  only  chanqe  beinq  that  the  vertical  distri- 
bution of  E  implicit  in  Fiq .  5(b)  is  sliqhtly  different  from  that  of 
the  constancy  of  E  everywhere  in  the  planetary  boundary  layer,  as  postu- 
lated by  Kaitala  (1974) „ 

E.   LARGE-SCALE  TURBULENT  HEAT  FLUX  ACROSS  THE  SEA-AIR  INTERFACE 

The  values  of  H„  represent  the  sensible  heatinq  rates  for  the  atmos- 
pheric layer  (8,10)  and  were  computed  for  each  qridpoint  alonq  the  four 
meridians.   The  values  of  Hr  in  qeneral  turned  out  to  be  larqe  and  posi- 
tive in  the  Southern  Hemisphere  reflectinq  the  computed  positive  values 

of  T  -  T   .   In  the  Northern  Hemisphere,  H_  was  predominantly  neqative 
g    x  r 

(cf.,  Figs.  7,..., 10)  and  therefore  contributes  negatively  to  the  net 
atmospheric  heat  budget  in  the  layer  (8,10). 

The  latent  heat  term  (E)  was  large  and  positive  relative  to  Hp  in 
the  Southern  Hemisphere  and  in  the  latitude  band  0-25N.   Regardless  of 
localized  variations  in  the  sign  of  Hr ,  combined  values  of  E  +  Hr  seemed 
to  give  reasonable  heating  rates  in  the  atmosphere  in  these  latitudes. 
On  the  other  hand,  the  combined  E  +  Hp  acts  as  heat- loss  contribution 
to  the  troposphere  in  Northern  Hemisphere  latitudes  <J>  >  30N  and  is  of 
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the  same  sign  as  the  sensible  heat-loss  contribution  (in  these  lati- 
tudes) applied  to  the  layer  (8,10). 

In  those  latitudes  (Southern  Hemisphere  and  the  Trade  Wind  Belt) 
where  E  +  H„  acts  as  a  positive  contribution  to  the  heat  budget  of  the 
tropospheric  column,  it  is  not  specified  in  this  study  precisely  where 
the  latent  heat  source  E  will  be  realized.   The  use  of  the  equations  of 
conservation  of  mass,  momentum  and  of  water- vapor  mass,  would  resolve 
the  latter  question,  but  this  was  considered  beyond  the  scope  of  the 
present  work. 

The  turbulent  flux  computations  of  Section  VI  for  H„  and  E  just  out- 
lined represent  minor  modifications  of  the  FNWC  turbulent  boundary 
layer  model  which  has  been  appended  to  the  major  concept  of  this  study 
namely  the  radiative  model.   The  purpose  for  inclusion  of  the  E,  Hp 
computations  has  been  primarily  for  consistency  tests,  that  is,  to 
determine  whether  net  radiative  flux  at  the  surface  is  of  the  proper 
magnitude,  relative  to  E  +  Hr . 


78 


VI.   MERIDIONAL  CROSS-SECTIONAL  DEPICTION 
OF  THE  HEAT-BALANCE  COMPUTATIONS 

A.  GENERAL 

The  general  design  of  this  section  is  to  utilize  all  of  the  computa- 
tional concepts  discussed  in  Sections  III,  IV,  and  V  in  the  computations 
for  a  single  time-step  in  the  heating  model  developed  for  use  in  the  FNWC 
prediction  model.   After  testing  FF2  and  ALB  with  the  corresponding 
values  of  Raschke  et  al  (1973) ,  it  was  decided  that  only  the  computations 
by  the  2/3-CL  parameterization  would  be  displayed  in  the  meridional 
cross-sections  (Figs.  6,7,8,9,10).   The  appropriate  calculations  were 
performed  at  each  gridpoint  of  the  four  meridians  used  for  presentation 
purposes  in  this  section.   However,  the  same  type  of  computations  using 
the  full-CL  parameterization  were  also  made,  but  not  presented  in  this 
section. 

B.  GEOGRAPHICAL   REPRESENTATION   OF    THE   HEAT-BALANCE   DISTRIBUTION 

The  FNWC  gridpoint  processed  analyses  for  0000GMT,  16  July  1974  were 
used  at  the  three  Pacific  cross-sections,  while  that  for  1200GMT,  16 
July  1974,  were  used  for  the  single  Atlantic  meridian.   This  was  done 
so  that  the  set  of  gridpoints  was  considered  to  be  subject  to  the  actual 
radiative- transfer  calculations  involved  for  these  specific  times. 
Figure  6  depicts  in  symbolic  language  the  key  to  the  computational  en- 
tries in  Figs.  7,8,9  and  10.   This  symbolic  list  presents  initially  the 
computations  made  at  each  gridpoint  (I, J)  radiative  sounding  data  in 
the  form  of  Table  I (b) .   The  computations  proceed  from  the  top  of  the 
troposphere  to  the  ocean  surface  and  include  the  sensible  and  latent 
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heat  transfers,  by  the  means  of  the  computational  procedure  of  Section  V. 
For  purposes  of  climatological  data  comparison,  Figs.  7,8,9  and  10  were 
developed  by  interpolating  gridpoint  results  to  whole  5°  latitude  incre- 
ments.  The  interpolation  routine  to  this  gridpoint  spacing  made  use  of 
the  Lagrangian  cubic  interpolation  scheme  (after  Spaeth,  1975) . 


(1-2)  (1-3)  (1-4)  (1-1)  (1-3)  (1-4) 

»*"*■'  «i      /i_o\  n_n  /i_/i\  «- 


•1    (1-2)  (1-3)  (1-4)         ^2     (2-1)  (2-3)  (2-4) 


(1-1)  (1-2)  (1-4)  (1-1)  (1-2)  (1-3) 

*!->      in  i  \    ii     o\    t->„A\  >i/ 


(6-1) 


•3    (3-1)  (3-2)  (3-4)         *4    (4-1)  (4-2)  (4-3) 


Finally  for  ease  in  reconciling  the  magnitudes  of  all  heat-transfer 
rates,  the  time-dependent  solar  disposition  rates  have  been  averaged  to 
24-hourly  rates. 

C.   EXPLANATION  OF  SYMBOLIC  TERMS 

lo   Cross-Section  at  Level  k=2  (Fig.  6) 

The  discussion  of  all  insolation  parameters  discussed  previously 
in  Section  IV  dealt  with  the  specific  time  of  day  that  corresponded  to 
the  hour  angle  h  for  the  instantaneous  time  t  under  consideration.   The 
incident  solar  insolation  dealt  with  is  then 

F(S)  =  S  (— )~2  Cos  z  (6-2) 

r 
m 

In  order  to  avoid  reference  to  specific  map  times  t,  the  instantaneous 
solar  hour-angles  were  h  =  5  5°,  10°,  3  5°,  and  35°,  respectively  for 
cross-sections  1,2,3  and  4,  as  depicted  in  Figs.  7,8,9  and  10. 

QAVE  represents  the  24-hour  average  of  F(2)  and  appears  as  the 
first  input  symbol  in  Fig.  6.   Its  value  is  considered  to  be  more 
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representative  climatologically  for  the  data  day  under  consideration 
than  F  ( 2 )  . 

QAVE  is  derived  by  the  formula 


QAVE  -  F(2)  £^_E  (6_3) 

Cos  z 


where 


Cos  z  =  [HSin(f)  Sin6  +  Coscf)  Cos6  SinH]/7T      (6-4) 


and 


H  =  arcCos  [-  tancj)  tan6]  (6-5) 


Here   6  =  21.5°   is  the  solar  declination  for  16  July  1974,  and  H  is 


the  hour  angle  at  local  sunset  at  latitude  (f) .   Cos  z  in  (6-3)  is  equal 
to  the  24-hour  average  cosine  of  the  zenith  angle,  Eq.  4-2.   The  24- 
hour  time  averaging  period  yielding  QAVE  gives  heating  results  consis- 
tent in  magnitude  with  the  terrestrial-flux  divergences,  which  change 
only  slightly  with  time  of  day.   This  conversion  to  expected  daily- 
averaged  solar  disposition  quantities  is  compatible  with  the  determina- 
tion of  a  heat  budget  for  the  given  date  (16  July  1974) . 
Other  parameters  needed  for  level  k=2  are 


,  ,  ,Cos  z,  ,,  ^v 

QREF  =  REF(t)  (Cqs   )  (6-6) 


where 

REF(t)  -  F(2)  -A26  -A610  -(IA10+IS10)        (6-7) 

REF(t)  is  the  instantaneous  solar  reflected  insolation  at  a  gridpoint 
and  QREF  is  its  24-hour  average,  assuming  that  the  planetary  albedo 
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remains  constant  for  the  24-hour  period,  although  it  was  calculated  for 
the  indicated  solar  time  t.   This  assumption  requires  that  the  cloud 
amounts  computed  at  the  indicated  synoptic  times  are  representative  of 
the  entire  day. 

The  same  principle  will  be  used  with  regard  to  all  other  solar 
parameters  in  the  conversion  from  time-dependent  values  at  solar  time  t 
to  24-hour  averaged  values.   Superior  bars  (   )  are  not  used  in  present- 
ing the  averaged  values  shown  in  the  cross-sectional  key,  Fig.  6.   The 
24-hour  average  earth- troposphere  system  balance,  BALT,  is  computed  from 

* 
BALT  =  QAVE  -  (QREF  +  F   )  (6-8) 

for  the  tropopause  level  k=2  at  the  indicated  latitude.   Net  terrestrial 

* 
fluxes,  such  as  F   ,  were  considered  to  be  constant  throughout  the  24- 
hour  period,  a  valid  assumption  if  the  cloud  cover  remains  quasi-con- 
stant for  the  period. 

2.   Cross-Section  in  Layer  (2,6) 

The  following  definitions  apply  for  the  layer  (2,6)  as  referred 
to  in  Fig.  6.   All  the  heat  transfers  shown  in  this  layer  are  assumed 
to  be  of  radiative  character  only,  as  in  the  corresponding  version  of 
the  FNWC  heating  model.   The  averaged  radiative  cooling  rate  is  given  by 

BAL  26  =  Q26  -  F26  (6-9) 


where  F26  =  terrestrial  cooling  rate  in 

layer  (2,6) 

and  Q26  =  daily  solar  absorption  in  layer  (2,6) 


is  defined  relative  to  A26(t)  by  a  cosine  transformation  similar  to 
Eq.  6-3. 
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3.  Cross-Section  in  Layer  (6,10) 

The  24-hour  average  radiative  cooling  in  the  layer  (6,10)  is 
given  by 

BAL  610  =  Q610  -  F610  (6-10) 

where  0610  is  the  sum  of  Q68  and  Q810 .    BAL  68  has  been  taken  as  one- 
half  of  BAL610.   However,  the  heat  balance  in  the  layer  (8,10)  has  been 
modified  for  inclusion  of  Hp,  as  given  by  Eq.  5-7. 

BAL810  =  Q810  -  F810  +  H   =  BAL68  +  H  (6-11) 

Here  Hr  is  the  sensible  heat-flux  convergence  in  the  layer  (8,10) . 
F610  is  the  sum  of  F68  and  F810,  the  terrestrial  flux  losses  for  the 
respective  layers. 

4.  Cross-Section  at  Air-Sea  Interface  (k=10) 

The  heat  balance  at  the  earth's  surface  (BALB)  is  as  defined 
in  the  following  equation 

BALB  =  QABG  -  F     -  (H   +  E)  (6-12) 

QABG  is  the  24-hour  average  insolation  absorbed  by  the  surface  defined 
by 


QABG  =  QABG(t)  (Cos  z/Cos  z)  (6-13) 

Hp  and  E  were  computed,  respectively,  from  Eqs .  5-7  and  5-13  and  con- 
sidered constant  for  the  computational  day,  16  July  1974. 
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D.   MERIDIONAL  CROSS-SECTIONS  OF  THE  VERTICAL  HEAT  BALANCE 
FOR  16  JULY  1974 

Figs.  7,8,9  and  10  as  previously  explained  represent  the  single  time 

step  of  heating  computations  for  each  of  the  four  meridians  used  in  this 

study.   The  four  figures  have  been  divided  into  (a)  tropical  results  and 

(b)  mid-to-high-latitude  results.   While  the  results  depicted  in  these 

cross-sections  are  exhibited  as  representing  daily-averaged  values,  they 

are  actually  based  upon  heat-computations  at  the  specific  map  times  of 

0000GMT  and  1200GMT  on  16  July  1974.   Therefore,  for  these  results  to 

be  meaningful  as  a  stepwise  part  of  the  heat  package  subroutine  of  FNWC, 

the  solar  radiative  absorption  and  reflectance  terms  would  have  to  be 

recoverable  as  a  function  of  GMT,  i.e., 


F(2,t)  =  CAVE   (Cos  z/Cos  z)  (6-14) 


REF(t)  =  OREF   (Cos  z/Cos  z)  (6-15) 

etc.   Thus  solar  disposition  terms  may  then  be  utilized  in  connection 
with  the  one-hour  time  stepwise  application  of  the  thermodynamic  equa- 
tion of  the  set  of  primitive  equations  used  in  the  FNWC  prediction 
process,  for  the  2/3-CL  parameterization  used  throughout  this  section. 

The  full-CL  heat-balance  cross-sections  were  also  computed,  but  have 
not  been  presented  here  for  the  sake  of  brevity.   However,  certain  key 
radiative  heating  (cooling)  rates  have  been  summarized  in  zonally- 
averaged  form  for  both  cloud  parameterizations  for  the  purposes  of  pre- 
senting the  mean  latitudinal  comparisons  (Section  VII)  of  the  purely 
radiative  heat  budget  of  the  FNWC  data  as  it  was  applicable  to  16  July 
1974. 
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VII.   THE  LATITUDINAL  DISTRIBUTION  OF  RADIATIONAL 
BALANCE  TERMS  OF  THE  OCEAN-ATMOSPHERE  SYSTEM 

A.   GENERAL 

The  latitudinally  distributed  results  of  the  purely  radiative  contri- 
butions, over  the  ocean-atmosphere  system,  are  presented  in  Figs.  12 (a, b) 
for  the  2/3-CL  and  13 (a, b)  for  the  full-CL  cases,  respectively.   These 
cross-sections  show  the  results  after  averaging  over  the  four  meridians 
considered  in  this  study.   The  results  are  displayed  in  the  format  of 
Fig.  11.   Radiative  transfers  Q68  and  Q810  have  been  combined  with  F68 
and  F810  of  Fig.  6  to  compute  layer-average  radiational  warming  rates  in 
the  layer  (6,10).   The  turbulent  transfer  heating  rates  E  +  H„  to  the 
atmosphere  (and  the  corresponding  cooling  rates  in  the  ocean)  which  were 
included  in  the  results  of  Figs.  7,8,9,  and  10  are  now  omitted  in  the 
radiative  composites  presented  in  Figs.  12  and  13. 

In  obtaining  latitudinally-distributed  means  of  radiative  heating 
rates  at  each  five-degree  multiple  of  latitude  (f) ,  all  values  of  each 
radiative  parameter  listed  in  Figs.  12  and  13  at  latitude  <f>  in  the  range 
20S,...,65N  (5  degree  increments)  were  simply  averaged  over  the  four 
meridians  providing  the  zonally-averaged  value  at  (J>.   At  (f>  =  65N  there 
was  only  one  contribution  to  the  latitudinal  average,  while  in  the  South- 
ern Hemisphere  latitudes  20,  15  and  10S,  only  two  values  (on  A  =  125W 
and  A  =  35W)  of  each  parameter  contributed  to  the  means  listed  in  Figs. 
12  and  13.   At  5S,  there  were  three  sets  of  radiative  parameters  enter- 


ing into  Q((J>).   Otherwise,  there  were  four  values  of  each  radiative 
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parameter  entering  into  the  computed  radiative  mean  values  Q((b)    of  the 
cross-sections.   Therefore,  near  the  northern  and  southern  boundaries 
of  Figs.  12  and  13,  the  listed  values  may  not  be  equally  representative, 
although  the  general  trend  should  be  reliable.   Fig.  14,  using  data 
extracted  from  Figs.  12  and  13,  tends  to  confirm  this  statement. 

Figs.  12,  13  and  14  are  also  useful  in  determining  the  relative 
merits  of  the  2/3-CL  model  as  compared  to  the  full-CL  model. 

B.   EARTH-TROPOSPHERE  SYSTEM  RADIATIONAL  BALANCE  SUMMARY 

The  latitudinally  averaged  distributions  of  the  three  key  radiative 

parameters,  R  ,  R  and  R  as  functions  of  latitude  for  both  the  2/3-CL 
s   a 

and  full-CL  cases  are  shown  in  Fig.  14.   The  radiative-balance  parameters 

were  defined  (after  Malkus,  1962)  as  R  ((})),  the  mean  radiative-energy 

transfer  rate  across  the  top  of  the  ocean-tropospheric  system  at  k=2 

(referred  to  as  BALT  in  Fig.  11);  R  ((b).    the  meridional  mean  radiative 

a 

cooling  rate  in  the  troposphere  (BAL26  +  BAL610  in  Fig.  11)  ;  and  R  ((b)  , 
the  mean  radiative  warming  (cooling)  rate  at  the  earth  surface  (BAL10 
in  Fig.  11).   The  relationship  between  these  radiative  parameters  is 


R  -  R  +  R  (7-1) 

s    a 


The  16  July  portion  of  Tables  X,  XI  and  XII  list  the  zonally-averaged 

values  of  R  ,  R  and  R  as  a  function  of  latitude  <b   for  the  2/3-CL  and 
s         a 
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the  full-CL  cases.   Weighted  averages  for  any  Q-value  are  then  computed 
using  the  following  cosine  weighting  scheme  (Eq.  7-2,  Spaeth,  1975) 


£  k.  *     a.. 

I  -t   <  Z  Ji  .  Cos  h 

0  =        1=1 3J=2 1 (7-2) 

vWt.  Avg.       18 

V    i 

/.  -7-  Cos  d) . 

1=1 


where  Q..  =  value  on  meridian  j  at  latitude  d> .    and  k.  is  the  number  k=l. 
31  1      x 

...,4  of  meridional  observations  available  for  the  meridional  averages 
at  (f> .  .   Note  that  i=l,...,18  corresponds  to  the  18  latitudes  (j) .  =  -20., 
,o.,65.  degrees.   Tables  X,  XI,  and  XII  also  depict  the  comparative  cross- 
seasonal  values  of  R  ,  R,  and  R  by  latitudes  for  both  2/3-CL  and  full-CL 

s  a 

parameterizations  for  16  April  1974,  16  July  1974,  and  16  October  1973, 

respectively. 

In  Fig.  14,  the  net-flux  terms  across  the  interfaces  k=2  and  k=10 

are  R  (<■))  and  R(d>) .   At  k=2  net-flux  values  obtained  from  the  2/3-CL 
s 

model  (denoted  R  , )  exceed  those  by  the  full-CL  model  (R  _)  at  all  lati- 
sl  s2 

tudes.   There  is  a  localized  maximum  value  of  both  R  ,  and  R  _  in  the 

si      s2 

latitude  zone  40  to  50N,  with  localized  minima  at  latitudes  45  and  55N. 

Table  IX  indicates  the  existence  of  maximum  CL  at  these  two  latitudes 

(45  and  55N)  suggestive  of  polar  frontal-zone  cloudiness  and  higher  than 

normal  global  albedo.   The  peak  hemispheric  values  of  R  ,  and  R  „  then 

si      s2 

occur  finally  at  latitudes  65N  in  this  study.   If  computational  data 
extended  further  poleward,  the  overall  maximum  R  could  possibly  occur 
at  latitudes  higher  than  65N,  in  view  of  the  high  declination  of  the 
sun  on  16  July  1974.   The  same  general  pattern  in  the  zonal  mean  radia- 
tive balance  R((f>)  at  sea-level  exists  for  both  the  2/3-CL  and  full-CL 
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Figure  14.   Radiational  balance  at  the  tropopause  (R  ) ,  at  the 
ocean  surface  (R) ,  and  in  the  tropospheric  column  (R  ) .   Solid 
lines  denote  computations  made  with  the  2/3-CL  (subscripted  "1") 
cloud  model  and  dashed  lines  correspond  to  the  full-CL  model, 
(subscripted  "2").   All  computations  for  for  16  July  1974  in  ly 
min-1  x  102. 
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cases.   However,  the  sea-level  maximum  in  the  40-50N  latitude  zone  is 
more  pronounced  than  that  at  latitude  65N.   It  is  to  be  noted  that  R 
for  the  2/3-CL  case  is  greater  than  R  over  the  entire  range  of  lati- 
tudes under  consideration  in  Fig.  14. 

In  particular  Fig.  14  for  16  July  1974,  shows  that  the  quasi-con- 
stant difference  R  ,  -  R  _  from  the  bottom  line  of  Table  X(b)  is 

si    s2 

.0346  ly  min   .   Likewise,  Fig.  14  shows  that  the  net-flux  difference 
R  (({>)  -  R_(<|>)  is  closely  approximated  by  .0303  ly  min   ,  the  result  of 
the  cosine-weighted  difference  from  Table  XI (b) . 

From  the  flux-divergence  equation  7-1  pertaining  to  the  radiative 
cooling  of  the  tropospheric  columns 


R   =  R_  -  R 
a    s 


the  values  of  R   and  R   turn  out  to  be  nearly  equal,  regardless  of 
a  J.      3.Z. 

cloud  configuration.   From  the  considerations  of  the  preceding  para- 
graph, the  cosine-weighted  mean  difference  R   -  R   is 

.  -1 


Ral  "  Ra2  =(Rsl  "  Rs2)-(R1  "  R2)=  -0043  ly  min 


which  also  results  from  Table  XII (b) .   This  is  a  relatively  small  differ- 
ence considering  the  mean  cooling  rates  involved.   Fig.  14  makes  this 

same  result  evident  upon  examination  of  the  curves  of  R  ,  and  R  n. 

al      a2 

There  are  consistent  features  in  both  the  R  ,  and  R  _  curves,  such  as 

al      a2 

the  minimum  values  in  the  latitude  zones  -10,  to  the  equator  and  20N  to 

30N,  which  make  the  general  results  on  R  ,  and  R  _  appear  credible. 

al      a2  cr 
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C.  CROSS-SEASONAL  EFFECTS 

Tables  X  and  XI  show  the  three  cross-seasonal  dates  affecting  R 

s 

and  R  centered  on  16  July  1974.   It  is  clear  that  R  ,  is  a  relative 

si 

cross-seasonal  minimum  at  all  (J)  <  5N  for  16  July  1974.   A  similar 
effect  applies  to  R   thrugh  the  same  time  sequence.   Intra-seasonal 
maximum  results  apply  poleward  of  ION  to  the  cross-seasonal  compari- 
sons of  R  at  the  earth's  surface  in  either  cloud-case. 

For  all  (J>  <_  5  degrees  of  latitude  in  Tables  X  and  XI,  there  is  a 
systematic  trend  toward  minimal  values  of  both  R   and  R,  at  16  July 

1974.   The  same  effect  applies  to  R  _  and  R^  on  this  date.   For  this 

s2      2 

limited  geographical  range  near  the  Southern  Hemisphere  border,  the 
16  October  date  having  a  solar  declination  6  =  -8.5°  latitude  corres- 
ponds to  the  largest  positive  values  of  both  R  and  R  in  either  cloud 
configuration.   These  cross-seasonal  results  appear  to  afford  quali- 
tative support  to  the  radiational  model  tested  here. 

D.  COMPARISONS  OF  NET  FLUX  AT  THE  TOP  OF  THE  MODEL  ATMOSPHERE 
WITH  SATELLITE  OBSERVATIONS 

Based  upon  observations  of  Raschke  et  al  (1973)  of  F2 (RAS)  and  of 

the  albedo  ALB (RAS)  over  the  latitude-range  of  this  study,  the  net 

flux  parameter  at  the  top  of  the  atmosphere  was  defined  for  the  NIMBUS 

III  data-period  (16-31  July  1969)  at  k=0 

RN(RAS)  =  QAVE*  T1^    (l-ALB(RAS))  -  F2(RAS)  (7-3) 

Raschke* s  computed  values  of  the  RN  for  this  data  period  used  a  slightly 
different  solar  constant  than  that  employed  in  this  study  (S  =  2.00  ly 
min-1) .    Hence  for  purposes  of  consistency  in  comparisons,  our  average 
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extraterrestrial  insolation,  QAVE* (2 .00/1. 92)  was  used  in  the  defini- 
tion of  RN(RAS)  in  Eq.  7-3. 

Values  of  RN(RAS)  were  computed  after  Raschke  at  each  five-degree 
latitude  intersection  on  the  four  meridians  where  comparisons  with  the 
present  model  computations  were  possible.   For  the  corresponding  radia- 
tive model  computations  of  this  study,  the  following  analog  to  Eq.  7-3 
was  adopted. 

RNMOD  =  QAVE*   *    (1-ALBMOD)  -  FF2  (7-4) 

FF2  represents  a  parameter  analogous  to  F2(RAS),  values  of  which 
were  compared  in  Section  III.Eo  on  a  latitudinal  basis.   Likewise, 
albedo  values  were  compared  on  a  latitudinal  basis  in  Section  IV. G. 

The  computational  values  of  RNMOD  for  16  July  1974  were  computed 
using  both  the  2/3-CL  and  the  full-CL  cloud  parameterizations .   After 
obtaining  averages  across  the  four  meridians,  the  latitudinal  compari- 
sons are  depicted  in  Table  XIII. 

Comparison  of  RN(RAS)  with  RNMOD  in  Table  XIII  shows  that  values 
of  RNMOD  are  somewhat  smaller  than  the  values  of  RN(RAS)  at  all  lati- 
tudes (J)  <  45N.   The  positive  difference 

RN(RAS)  -  RNMOD       <f)  <  45N 

occurs  because  of  the  consistently  small  albedos  of  the  Raschke 
observations.   This  difference  is  considerably  larger  for  the  full-CL 
case,  when  the  model  albedos  were  excessively  high. 


105 


R  -values  for  the  earth-troposphere  have  been  extracted  from  Table 
s 

X  and  appear  in  columns  four  and  five  of  Table  XIII.   Recall  that 


R   =  pAVE  (1-ALBMOD)  -  F  *  (7-5) 

s    '  2 


In  comparing  Eqs.  7-4  and  7-5,  it  is  clear  that  the  incoming  insolation 
used  in  the  computation  of  R  has  been  reduced  by  4%  of  an  amount  which 


is  approximately 


-1 


QAVE  -  QREF  =  .4148  ly  min  (7-6) 


when  cosine-averaged  over  all  latitudes  considered  in  this  study, 
Table  XIII  indicates  that 

.  -1 


R  ,  -  RN (MOD- 2/3)  =  .025  ly  min 
si 


holds  for  the  2/3-CL  case.   This  difference  may  be  reconciled  with  the 
net-flux  differences  of  Eqs.  7-5  and  7-4,  through  the  differenced 
equation 


R    -  RN(MOD-2/3)  =  (FF2  -  F  *)  -  . 04 (QAVE  -  QREF)  (7-7) 


which  requires  that 


(FF2  -  F  *)  =  .025  +  .04 (QAVE  -  QREF)  (7-8) 

in  ly  min      for  the  2/3-CL  case. 

The  result  of  Eq.  7-8  shows  the  average  I R- flux  divergence  across 
the  layer  k=2  to  k=0  by  the  F   model  of  Eq.  3-5.   The  water-vapor  and 
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CO   absorber  masses  of  the  stratosphere  are  therefore  not  quite  in  radia- 
tive equilibrium,  because  these  same  absorber  masses  do  not  have  a  com- 
pensating solar  absorption  effect  in  the  present  model.   The  4%  solar 
absorption  by  ozone  and  oxygen  are  not  sufficient  to  compensate  for  the 
added  water  vapor-CO   flux-divergence  in  the  terrestrial  waves. 
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LAT. 
2  OS 
15 
10 

5 

0 

5 
10 
15 
20 
25 
30 
35 
40 
45 
50 
55 
60 
65 

Wt.Avg. 


Mid- 

-Seasonal 

Dates 

(a) 

(b) 

(c) 

16  April  1974 

16  July 

1974 

16  October  1973 

R  , 
si 

RS2 

R  i 
si 

s2 

Rsl 

s2 

0229 

-.0469 

-.1012 

-.1124 

.0764 

.0384 

0119 

-.0418 

-.0817 

-.1010 

.0785 

.0400 

0002 

-.0358 

-.0546 

-.0772 

.0803 

.0430 

0171 

-.0235 

-.0468 

-.0656 

.0814 

.0480 

0353 

-.0079 

.0138 

-.0154 

.0893 

,0567 

0283 

.0243 

.0297 

-.0062 

.0696 

.0388 

0198 

-.0364 

.0521 

.0145 

.0484 

.0164 

0225 

-.0335 

.0635 

.0211 

.0170 

-.0030 

0204 

-.0368 

.0627 

.0136 

-.0011 

-.0322 

0256 

- .0126 

.0825 

.0380 

-.0281 

-.0588 

0547 

.0293 

.0945 

.0560 

-.0573 

-.0840 

0488 

.0301 

.1018 

.0607 

-.0743 

-.0916 

0540 

.0296 

.1254 

.0888 

-.0916 

-.1036 

0427 

.0096 

.0892 

.0612 

-.1131 

-.1212 

0055 

-.0362 

.1461 

.1038 

-.1347 

-.1395 

0016 

-.0146 

.1063 

.0766 

-.1307 

-.1302 

0428 

-.0591 

.1603 

.1315 

-.1525 

-.1493 

0089 

-.0374 

ol990 

.1640 

-.1169 

-.1096 

0027 

-„0153 

.0536 

.0190 

-.0081 

-o0320 

Table  X.   Cross-seasonal  comparison  of  zonally-averaged  values  of  R 
as  a  function  of  latitude.   The  averaged  radiational  values  for  the" 
2/3-CL  cloud  parameterization  are  denoted  with  a  subscripted  1  and 
the  full-CL  values  by  a  subscripted  2.   (All  values  in  ly  min~l) . 
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Mid- 

Seasonal  Dates 

(a) 

(b) 

(c) 

16  April 

1974 

16  July 

1974 

16  October  1974 

LAT. 

Ri 

R2 

Ri 

R2 

Ri 

R2 

2  OS 

.0813 

.0514 

.0598 

.0547 

.1742 

.1329 

15 

.1015 

.0681 

.0954 

.0852 

.1845 

.1438 

10 

c1299 

.0957 

.1277 

.1134 

.1911 

01524 

5 

.1438 

.1049 

.1397 

.1255 

.1756 

.1352 

0 

.1164 

.1376 

.1741 

.1508 

.2117 

.1787 

5 

.1858 

.1294 

.1773 

.1464 

.2017 

.1722 

10 

.1911 

.1460 

.2042 

.1725 

.1886 

.1598 

15 

.1933 

.1470 

.2216 

.1866 

.1657 

.1465 

20 

.2018 

.1569 

.2373 

.1987 

.1611 

.1365 

25 

.2079 

.1769 

.2580 

.2225 

.1490 

.1282 

30 

.2153 

.1917 

.2016 

.2286 

.1179 

.0988 

35 

.1889 

.1684 

.2774 

.2429 

.0870 

.0712 

40 

.1824 

.1549 

.2920 

.2598 

.0791 

.0711 

45 

.1670 

01343 

.2474 

.2224 

.0547 

.0507 

50 

.1256 

.0952 

.2751 

.2313 

.0341 

.0342 

55 

.1062 

.0850 

.1824 

.1447 

.0112 

.0151 

60 

.0818 

.0603 

.2724 

.2359 

-.0040 

.0017 

65 

.0614 

.0206 

.2637 

.2154 

.0060 

.0119 

Wt.  Avg. 

.1684 

.1333 

.2119 

.1816 

.1370 

.1155 

Table  XI.   Cross-seasonal  comparison  of  zonally-averaged  values  of 
R  as  a  function  of  latitude.   The  averaged  values  for  the  2/3-CL 
case  are  denoted  with  a  subscripted  1  and  the  full-CL  case  values 
by  a  subscripted  2.   (All  values  in  ly  rain--1-). 
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Mid- 

-Seasonal  1 

Dates 

(a) 

(b) 

(c) 

16  April 

1974 

16  July 

1974 

16  October  1973 

LAT. 

R  i 
al 

a2 

al 

a2 

R  i 
al 

a2 

2  OS 

-.1042 

-.0983 

-.1613 

-.1671 

-.0900 

-.0942 

15 

-.1134 

-.1099 

-.1772 

-.1863 

-.1061 

-.1039 

10 

-.1297 

-.1315 

-.1824 

-.1906 

-.1107 

-.1093 

5 

-.1267 

-.1284 

t.1869 

-.1912 

-.0941 

-.0872 

0 

-.0811 

-.1455 

-.1604 

-.1663 

-.1225 

-.1219 

5 

-.1575 

-.1671 

-.1476 

-.1527 

-.1322 

-.1334 

10 

-.1713 

-.1824 

-.1522 

-.1582 

-.1402 

v. 1434 

15 

-.1708 

-.1805 

-.1582 

-.1657 

-.1488 

-.1495 

20 

-.1814 

-.1937 

-.1749 

-.1853 

-.1621 

-.1687 

25 

-.1823 

-.1895 

-.1756 

-.1846 

-.1772 

-.1869 

30 

-.1606 

-.1622 

-.1673 

-.1727 

-.1753 

-.1828 

35 

-.1401 

-.1383 

-.1757 

-.1824 

-.1613 

-.1628 

40 

-.1284 

-.1253 

-.1667 

-.1711 

-.1707 

-.1748 

45 

-.1243 

-.1247 

-.1595 

-.1612 

-.1679 

-.1718 

50 

-.1311 

-.1314 

-.1291 

-.1277 

-.1688 

-.1736 

55 

-.1046 

-.0996 

-.0762 

-.0684 

-.1419 

-.1453 

60 

-.1246 

-.1194 

-.1123 

-.1046 

-.1486 

-.1511 

65 

-.0703 

-.0580 

-.0648 

-.0517 

-.1226 

-.1215 

Wt„  Avg„ 

-.1457 

-.1486 

-.1585 

-.1628 

-.1450 

-.1474 

Table  XII.   Cross-seasonal  comparison  of  the  zonally-averaged  values 
of  R  as  a  function  of  latitude.   The  averaged  values  for  the  2/3-CL 
are  denoted  with  a  subscripted  1  and  the  full-CL  cases  values  by  a 
subscripted  2.   (All  values  in  ly  min-1) . 
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Ocean 

RNMOD 

Troposph 

ere 

LAT. 

2/3-CL 

Full-CL 

R  , 
si 

RS2 

RN(RAS, 

2  OS 

-.1144 

-.1256 

-.1012 

-.1124 

-.0594 

15 

-.0928 

-.1121 

-.0817 

-.1010 

-.0362 

10 

-.0618 

-.0844 

-.0546 

-.0772 

-.0038 

5 

-.0554 

.0287 

-.0468 

-.0656 

.0182 

0 

-.0042 

-.0345 

.0138 

-.0154 

.0616 

5 

.0108 

-.0252 

.0297 

-.0062 

.1194 

10 

.0374 

-.0002 

.0521 

.0145 

.1372 

15 

.0572 

-.0102 

.0635 

.0211 

.1430 

20 

.0577 

.0086 

.0627 

.0136 

.1500 

25 

.0632 

.0348 

.0825 

.0380 

.1426 

30 

.0873 

.0488 

.0945 

.0560 

.1484 

35 

.0958 

.0546 

.1018 

.0607 

.1494 

40 

.1148 

.0782 

.1254 

.0888 

.1561 

45 

.0777 

.0504 

.0892 

.0612 

.1164 

50 

.1244 

.0819 

.1461 

.1038 

.0965 

55 

.0617 

.0317 

.1064 

.0766 

.0997 

60 

.1323 

.1034 

.1603 

.1315 

.0715 

65 

.1436 

.1082 

.1990 

.1640 

.0708 

Wt.  Avg. 

.0288 

.0027 

.0536 

.0190 

.0858 

Table  XIII.   Comparison  of  the  ocean-troposphere  net  radiation  R  at 
k=2,  with  the  model  computed  net  radiation  at  the  top  of  the 
atmosphere  RNMOD  and  with  climatology  net  radiation  at  the  top  of 
the  atmosphere  from  Raschke  et  al  (1973)  RN(  RAS)  .  R  and  RNMOD  are 
shown  for  both  cloud  cases. 
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VIII.   ZONALLY-AVERAGED  TROPOSPHERIC  AND 

OCEANIC  HEAT  BUDGETS  FOR  16  APRIL  1974 

A.   THE  TROPOSPHERIC  HEAT  BUDGET 

By  averaging  over  the  tropospheric  columns  for  the  four  meridians 
displayed  in  Figs.  7,  8,  9  and  10,  latitude  by  latitude,  the  zonal 
distribution  of  E  +  Hp  results  were  obtained  and  are  presented  in 
Table  XIV.   Recall  that  R  ,  the  tropospheric  radiative  net  cooling  is 
listed  in  Table  XII.   The  resulting  heating  rate  of  the  troposphere 
may  then  be  expressed  as  a  function  of  <J)  by  the  right  side  of  Eq.  8-1 

Q    +  S   =  R   +  (E  +  Hr)  (8-1) 

va    a    a         T 

LAT.  20S     15      10       5  0  5      10      15      20 

E+H  .7044   .3784   .2014   .0273  .0595  .0592   O0669   .1296   .0761 

LAT.  25      30      35      40  45  50      55      60      65 

E+Hr  .0852  -.0066  -.0180  -.0411  -.0047  -.0344  -.0358  -.0092  -.0303 

Table  XIV.   Zonally-averaged  values  of  E  +  Hr  as  a  function  of  latitude, 
(All  values  in  ly  irtin   )  . 


S   is  the  storage  heating  rate  of  the  column,  and  Q    is  the  re- 
a  *va 

quired  flux-divergence  of  heat  compatible  with  the  heat  balance  equa- 
tion for  the  column,  for  the  observed  values  of  R  ,  E  +  Hn  and  S  . 

a       l       a 

Figure  15  depicts  the  three  heating  rates  (R  ,  (E  +  Hp)  and  (Q   +  S  ) ) 

a        1         va    a 

as  a  function  of  latitude  for  the  two  cloud-model  cases. 
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20S 
\ 


(E  +  Hr) 


(Qva 


(Ova  +  Sa), 


20S 


0 


10  20 

LATITUDE 


60  65N 


Figure  15.   Latitudinal  distribution  of  R  ,  E  +  Hp ,  and  of  Q   +  S 
as  a  function  of  latitude  for  the  two 

cloud-model  cases.   Values  of  E  +  Hp  at  15  and  20S  were  off  the 
upper  scale  of  the  ordinate  and  may  be  found  in  Table  XIV.   (Solid 
curves  denote  2/3-CL  case,  dashed  curve  full-CL) . 
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As  noted  in  Section  VII. B.,  the  difference  R  ,  -  R  ^  is  almost  negli- 

al    a2 

gible  in  this  study,  and  there  is  little  variation  as  a  function  of 

latitude.   The  distribution  of  E  +  H„  as  a  function  of  latitude  is  also 

independent  of  cloud-model.   Therefore,  Fig.  15  shows  only  one  curve  for 

E  +  H„.   The  curve  of  (Q   +  S  ) ,  is  shown  in  Fig.  15,  as  the  residual 
T  va    a  1 

of  the  right  side  of  Eq.  8-1  for  the  2/3-CL  case. 

For  ease  of  interpretation  in  Fig.  15,  the  curve  of  (Q   +  S  )   for 

va    a  2 

the  full-CL  case  was  not  plotted.   This  curve  follows  closely  parallel  to 

that  of  (Q   +  S  ) .,  being  slightly  less  than  (Q   +  S  )  _  in  the  latitude 
xva    a  1  va    a  1 

range  20S  to  45N.   From  45N  to  65N  the  curve  of  (Q   +  S  )„  exceeds  the 

va    a  2 

curve  of  (Q   +  S  ) ,  by  a  small  amount.   This  crossover  effect  is  to  be 
va    a  1 

expected  by  examination  of  the  corresponding  crossover  in  the  R   and 
R   curves  at  these  higher  latitudes.   As  is  seen  in  Fig.  15,  the  net 

3.2. 

heating  curve  Q   +  S   follows  the  same  basic  pattern  as  the  E  +  H„ 
*va    a  T 


curve,  since  R  is  relatively  constant  with  latitude, 
a 


-1 


The  cosine-weighted  value  0   +  S   =  -.0512  ly  min  "  for  the  average 

"va    a 


tropospheric  column.   Q   has  not  been  computed  in  this  study;  however, 


if  Q   +  S   is  attributed  solely  to  mean  storage-cooling  S   of  the  tropo- 
va    a  a 

sphere,  the  daily  storage  rate  corresponds  to  a  temperature-change  rate 
given  by 


—         R  +  (E  +  H  ) 

lk]    =   4a    APmb X  144°  °C/day        (8"2) 


with  AP  mb  =  800  mb  in  the  troposphere.  The  resultant  cooling  rate  by 
(8-2)  in  the  troposphere,  considering  that  zero  lateral  flux-divergence 
applies,  is  therefore  approximately 

(21)  =  _.59°C  (day)"1 
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2 

averaged  over  the  mean  cm  tropospheric  column. 


B.   THE  LATITUDINALLY-AVERAGED  HEAT  BUDGET  OF  THE  OCEAN 

The  latitudinally-averaged  heat-budget  of  the  ocean  water-mass 
obeys  the  following  thermodynamic  relationship  for  both  the  2/3-CL 
and  full-CL  cases. 


Q   +  S   =  R  +  (-(E  +  Hr))  (8-3) 

vo     o  I 


Note  that  Eq.  8-3  is  the  analog  of  Eq.  8-1  for  the  ocean-mass,  with 
- (E  +  H„)  representing  the  average  turbulent  cooling  rate  for  the 
water-mass  by  mixing  across  the  sea-air  interface. 

The  two  terms  of  the  right  side  of  Eq.  8-3  have  been  computed  for 
each  meridian  and  then  averaged  zonally,  with  R  and  E  +  H„  shown  in- 
Table  XI  and  XIV,  respectively.   0   +  S   has  been  computed  as  the 

xvo    o  c 

residual  of  these  two  terms  at  each  latitude  and  is  graphed  in  Fig.  16 

for  both  cloud  cases. 

The  primary  ocean-mass  heating-function  is  R,  the  radiative  heating 

rate  at  the  ocean-air  interface.   As  discussed  in  Section  VII,  R  >  R 

for  all  latitudes  considered  in  this  study.   Therefore,  according  to 

Eq.  8-3,  (Q   +  S  ),  must  be  greater  than  (Q   +  S  ) _  by  the  same 
vo    o  1  vo    o  2 

amounts  since  - (E  +  Hp)  is  constant  for  both  cloud  cases.   This  means 
that  the  internal  heating  rate  of  the  ocean-mass,  as  given  in  Eq.  8-3, 
is  greater  in  the  2/3-CL  case  than  in  the  full-CL  case  for  all  latitudes 
considered  in  this  study.   Q    is  the  required  mean  oceanic  heat-flux 
divergence  at  latitude  (f>,  and  S   is  the  heat-storage  in  the  oceanic 
column. 
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Figure  16.   Latitudinal  distribution  of  R,  - (E+Hn) ,  and  of  (Q  +S  ) 
.  1  vo  o 

as  a  function  of  latitude  for  the  two  cloud-model  cases.   The 

value  of  -(E+Hr)  at  15  and  20S  were  off  the  lower  scale  of  the 

ordinate  and  may  be  inferred  in  Table  XIV.   (Solid  curves  denote 

2/3-CL  case;  dashed  or  dot-dashed  curves  denote  full-CL  case) . 
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The  cosine-weighted  average  of  (Q   +  S  ) ,  over  the  latitude  range 

vo     o  1 

of  Fig.  16  is 

Jo        +  S  ),    =  .0949  ly  min" 
vo     o  1 

(cf.,  Tables  XI,  XIV)  which  corresponds  to  a  mean  heating  rate  in  the 
water-mass  column.  It  should  be  noted  that  this  heating  rate  is  con- 
siderably larger  (and  of  opposite  sign)  than  the  corresponding  atmos- 
pheric effect. 
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IX.   CONCLUSIONS 

The  present  work  continues  the  use  of  the  radiational  two-layer 
cloud  model  adapted  from  Spaeth's  mid-winter  study  to  FNWC  gridpoint 
data  valid  over  the  four  meridians  for  16  July  1974.   The  use  of  the 
2/3-CL  parameterization,  as  compared  to  the  full-CL  after  Smagorinsky 
(1960) ,  for  the  specification  of  cloud-layer  amounts  gave  the  best 
verification  for  the  computed  global  albedo  as  a  function  of  latitude 
when  compared  with  observations  from  NIMBUS  III  satellite  climatology 
comparable  to  that  of  16  July  1974.   The  use  of  the  2/3-CL  parameteri- 
zation also  gave  good  verification  of  computed  long-wave  flux  to 
space  when  compared  to  the  latitudinal  distribution  of  the  same  element 
as  given  by  satellite  climatology. 

The  chief  shortcoming  of  the  computed  global  albedo  versus  the 
satellite  measured  albedo  occurred  in  tropical  latitudes  (20S  -  20N) . 
Here,  it  appears  that  cloud-elements  do  not  behave  as  reflecting  sheets 
for  insolation  but  rather  as  focusing  devices  through  the  intra-cloud 
spaces.   The  use  of  50-60%  model  cloud  reflectances  seems  excessive  in 
these  latitudes.   As  a  practical  alternative  it  may  be  reasonable  to 
tune  the  cloud  reflectances  until  agreement  of  the  computed  global 
albedos  with  satellite  measurements  is  achieved  in  the  latitude  range 
20S  -  20N. 

In  the  topics,  all  analyses  are  subject  to  data-sparsity  and  hence 
are  largely  dependent  upon  climatology  and  a  relatively  few  radiosonde 
reporting  stations.   An  apparent  weakness  in  this  study  is  that  of 
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overspecifying  cloud-amounts  through  the  use  of  Smagorinsky  type 
formulas.   This  could  be  a  result  of  inadequately  resolved  relative 
humidities,  which  in  general  might  well  be  lower  than  inferred  by 
the  FNWC  analyses. 

One  of  the  major  positive  results  of  the  radiative  model  was 
the  simulation  of  the  proper  sense  of  the  cross-seasonal  radiative 
heat  balances  at  the  top  of  the  earth-atmosphere  system.   A  similar 
positive  result  applied  to  cross-seasonal  radiative  balance  at  the 
ocean  surface  centered  on  16  July  1974. 
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